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Abstract 



In this thesis, wc consider the dynamics of vortices in the easy plane insulating 
ferromagnet in two dimensions. In addition to the quasiparticle excitations, 
here spin waves or magnons, this magnetic system admits a family of vortex 
solutions carrying two topological invariants, the winding number or vorticity, 
and the polarization. 

A vortex is approximately described as a particle moving about the system, 
endowed with an effective mass and acted upon by a variety of forces. Clas- 
sically, the vortex has an inter-vortex potential energy giving a Coulomb-like 
force (attractive or repulsive depending on the relative vortex vorticity), and 
a gyrotropic force, behaving as a self-induced Lorentz force, whose direction 
depends on both topological indices. 

Expanding semiclassically about a many-vortex solution, the vortices are quan- 
tized by considering the scattered magnon states, giving a zero point energy 
correction and a many-vortex mass tensor. The vortices cannot be described 
as independent particles — that is, there are off-diagonal mass terms, such as 
^MijViVj, that are non-negligible. 

This thesis examines the full vortex dynamics in further detail by evaluating 
the Feynman- Vernon influence functional, which describes the evolution of the 
vortex density matrix after the magnon modes have been traced out. In addition 
to the set of forces already known, we find new damping forces acting both 
longitudinally and transversely to the vortex motion. The vortex motion within 
a collective cannot be entirely separated: there are damping forces acting on one 
vortex due to the motion of another. The effective damping forces have memory 
effects: they depend not only on the current motion of the vortex collection but 
also on the motion history. 
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Chapter 1 



Introduction 



In a wide variety of systems, there exist vortices, high energy states nonetheless 
significant in system dynamics at low temperatures. Despite its high energy, a 
vortex can nonetheless form via tunneling processes or at a boundary with only 
a small energy barrier. They are exceptionally stable, arguable topologically, 
and, in fact, can only be destroyed if one meets its 'anti-vortex' or, equivalently, 
annihilates at a boundary (where it has met its image vortex) . Cooling a system 
down vortex-free is non-trivial, and, in general, we retain a low density of vortex 
states down to the lowest temperatures. 

Quantum vortices were first proposed in the 1950's in superfluid helium to ex- 
plain the decay of persistent currents. Since then, they have been proposed and 
measured in, for example, superconductors and a variety of magnetic systems. 
The dynamics are well described phenomenologically as a point-like particle in 
2D (or as a line in 3D) endowed with an effective mass and acted upon by a 
variety of forces. Microscopic derivations of the particle properties of a quantum 
vortex have been plagued by decades of debate and controversy. A recent resur- 
gence in debates began in the 1990's concerning the so-called Magnus force, a 
force borrowed from classical fluids acting perpendicular to the velocity. Ao and 
Thouless^ claimed that in superfluid helium (He II) there is a universal form 
of this force, independent of quasiparticle scattering. Others argue that there 
should be, in addition to the bare Magnus force, a tranverse damping force, 
reinforcing or opposing the Magnus forc e 

In this thesis, we consider a relatively simple magnetic system, a 2D insulat- 
ing ferromagnet with easy plane anisotropy, admitting a family of topologically 
stable vortices. We derive microscopically the vortex effective mass and, in 
addition to the previously reported gyrotropic force, the magnetic analogue to 
the Magnus force, and inter-vortex Coulomb-like forces, we derive a variety of 
vortex damping forces. We find both the usual longitudinal damping force and 
a transverse damping that acts in combination with the gyrotropic force. A 
transverse damping force has not yet been considered in a magnetic system. In 
fact, all treatments of the dissipative motion of a vortex have been phenomeno- 
logical, with the exception of Slonczewski'sS treatment with which we compare 
results in Chapter 2| A collection of vortices cannot be considered as a set of 
independent particles — they have mixed inertial terms and damping force terms. 
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Chapter 1. Introduction 



We first review a few symraetry arguments for the existence and stability of 
vortex solutions. Besides revealing the similarity between vortices from vari- 
ous systems, we find that vortices are an example of a more general family of 
topological solitons. 

We then briefly discuss the early work on quantizing solitons by the relativistic 
field theorists, focussing rather on the techniques than the various specific con- 
tributions. Note that we will use many of these techniques for quantizing the 
vortex in the easy plane magnetic system. 

Next, we discuss briefly solitons in condensed matter systems and the excit- 
ing new phenomena found there. For example, by examining the conducting 
polymers, fractional charge was first predicted and observed. 

Returning specifically to vortices, we briefly discuss the controversy in the mi- 
croscopic derivation of the equations of motion for a superfluid vortex. This 
will introduce the variety of forces wc should expect to act on a collection of 
vortices. Switching to magnetic systems, we find that despite the ease of direct 
experimental observation and simplicity of calculations not much work has been 
done here. 

Finally, we introduce in detail the magnetic system under consideration. The 
symmetry of the system admits topologically stable vortices and gapless quasi- 
particles. The purpose of this thesis is to separate the quantum dynamics of 
the vortices from the effects of the perturbative quasiparticles, here magnons. 



1.1 Symmetry breaking 

Symmetry plays a crucial role in science and we strive to discover and exploit 
the symmetries of the laws of nature (Galilean or Lorentz invariance, gauge 
invariance, etc.). However, we find that the symmetry of physical states may 
be a smaller subset of the full symmetry in which it resides. For example, in a 
Heisenberg ferromagnet, we find a system of spins free to lie in any direction in 
3D, preferring to align parallel to one another, however, in the absence of any 
magnetic fields, with no preference of which direction along which to lie. The 
ground state then chooses at random along what direction to align. 

A system with a degenerate ground state is forced to spontaneously choose one 
state amid the degeneracy, an example of spontaneously broken symmetry. A 
discrete degeneracy is found in the problem of a field residing in a double well 
potential (as in Figure ITtI left), or, more generally, an n-well potential. A 
continuous degeneracy in a system has a continuum of minima in the potential 
(as, for example, in Figure [TTTl right). The ferromagnet is an example of a 
system with a continuum of ground states, except that here, the potential is 
completely fiat: there is no preference at all between directions. 



1.1. Symmetry breaking 
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Figure 1.1: Left, an example potential of a ID field with a doubly degener- 
ate ground state; right, an example potential of a 2D field with a 
continuum degeneracy in its ground state. 



In general, different regions of a sample may choose different degenerate states 
or may even lie in an excited state. A mapping of the state taken across the 
sample, in all its available degrees of freedom, is called the order parameter. In 
a Heisenberg spin system, this is simply the spin vector in 3D as a function of 
position in the sample. The order parameter here can be mapped onto a unit 
sphere — a path along the sample is then traced as a path on the surface of the 
sphere. For a spin system confined to lie in the plane, the so-called XY model, 
the order parameter is mapped onto the unit circle. 

Incidentally, the order parameter in superfluid helium II can also be mapped 
onto the unit circle so that it is topologically equivalent to the XY model. This 
does not mean, however, that the dynamics of the vortices in each system should 
be the same, but, rather, only that the topology of vortices is identical in the 
two systems. 

If a system possesses discrete symmetries, to pass from one ground state to an- 
other there must be some transition region, or domain wall, separating different 
states. This domain wall, sometimes called a kink, is an example of a quasi-lD 
soliton. 

For a continuous symmetry, we can imagine similar cases where certain regions 
are forced out of a ground state. As a simple example, consider the XY spin 
model. If the spins choose to nearly align along the boundary, turning very 
slowly so as to always radiate outward, as we near some central region the 
spins are less and less ferromagnetically aligned and, further, there is a point 
discontinuity at the very center (see Figure [L^ . 

If we follow a path surrounding the vortex in order parameter space, that is along 
the unit circle, we find we must wrap around the unit circle once. This vortex 
is called a topological soliton with single wrapping number or vorticity. In this 
example, no matter how we smoothly deform the spins, we cannot continuously 
deform away this wrapping of the unit circle. We say that it is homotopically 
distinct from a zero winding path, or more simply a point. 
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Figure 1.2: A 2D XY-ferromagnet with a vortex connecting the degeneracy of 
spin directions. The central red dot signifies the point of disconti- 
nuity. 




Figure 1.3: A magnetic vortex formed by Heisenberg spins can be continuously 
deformed away by expanding about a patch of the unit sphere not 
covered by the vortex path, shrinking the vortex to a point. 

There are vortices with higher winding numbers, always integral to ensure con- 
tinuity. Each family of solutions corresponding to a certain winding number is 
topologically stable. That is, there exists no homotopy, or continuous mapping, 
between solutions of differing winding numbers. 

There are systems that admit vortices for which this topological stability is not 
guaranteed, and are thus not called topological solitons. Consider a general 
vortex residing on a sample for which the order parameter maps onto a unit 
sphere (Figure [l~5|l . The vortex is homotopically equivalent to a point (that is, 
a region with constant ground state) since we can imagine continuously shrinking 
the vortex away. In real space, this is equivalent to the ability of the spins to 
unwind, that is, all the spins twisting to all lie parallel to one another. Note that 
this unwinding is a special feature of the isotropy of the system. Although such 
a soliton does not possess topological stability, the entire plane must unwind, a 
macroscopic number of spins in the magnetic vortex case, so that the soliton is 
still essentially stable. 

The vortices considered in this thesis have an order parameter lying on the unit 
sphere, however, with a higher potential at the north and south poles. They 
are very similar to the XY vortex shown in figure 11.21 except that the spins 
are not entirely restricted to lie in the plane and, at some energy expense to 
restore continuity, the spins twist out of plane at the vortex center choosing 
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Figure 1.4: A vortex with +1 winding in a 2D Heisenberg ferromagnet with 
spins lying prcferentiaUy in the plane. 




spontaneously between the two possible perpendicular directions in which to 
twist. This direction is a second topological invariant of the vortices and is 
termed the polarization. An example of a vortex with unit winding number, or 
vorticity, and polarization out of the page is shown in Figure 11.41 There exist 
also zero polarization vortices lying entirely in the plane. 



1.2 Classical Solitons 



We found that vortices are examples of a topological solitons. Generally, a 
soliton is a finite energy localized solution of a wave equation, satisfying strict 
stability conditions under collisions with other soliton solutions 

tSee, for instance, the excellent book by Rajaraman^ on the quantization of solitons for 
a rigorous definition of a soliton. 
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The first reported soliton was in 1834 by John Scott Russell ^"^ in the Union 
Canal near Edinburgh (see Figure 



I was observing the motion of a boat which was rapidly drawn along a 
narrow channel by a pair of horses, when the boat suddenly stopped - 
not so the mass of water in the channel which it had put in motion; 
it accumulated round the prow of the vessel in a state of violent 
agitation, then suddenly leaving it behind, rolled forward with great 
velocity, assuming the form of a large solitary elevation, a rounded, 
smooth and well-defined heap of water, which continued its course 
along the channel apparently without change of form or diminution 
of speed. I followed it on horseback, and overtook it still rolling on 
at a rate of some eight or nine miles an hour, preserving its original 
figure some thirty feet long and a foot to a foot and a half in height. 
Its height gradually diminished, and after a chase of one or two miles 
I lost it in the windings of the channel. Such, in the month of August 
1834, was my first chance interview with that singular and beautiful 
phenomenon which I have called the Wave of Translation. 



He went on to build a 30' wave tank in his back garden in which to conduct 
further experiments on his "waves of translation" . 

In physics, there are the familiar optical solitons, with which demonstrations 
of long haul, low bit-error-rate transmissions have been made. In optics, a 
soliton is a localized EM wave with much higher power than a traditional optical 
signal. However, as opposed to regular low power optical transmissions, an 
optical soliton does not suffer dispersion, so that a signal is not distorted when 
transmitted over large distances. 

A soliton is usually a solution to a partial differential equation in which com- 
peting non-linear terms cooperate to create a self-reinforcing large amplitude 
solution. For instance, for a non-linear dissipative system, ordinarily, wave so- 
lutions are dispersive, that is, different k modes separate, and dissipative, energy 
spreads in real space. For these special soliton solutions the two mechanisms 
can act in opposition so that the net result is a non-dispersive, non-dissipative 
wave. 

More specifically, however, a vortex is an example of a topological soliton. These 
exist, not because of finely balanced non-linear terms in the equations of motion, 
but rather due to a degenerate freedom in the boundary conditions entailing the 
existence of homotopically distinct solutions (that is, solutions for which there 
is no continuous deformation from one to another). 
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Figure 1.6: An illustrative potential of a one dimensional particle. A soliton is 
analogous to the second minimum at a; = c. 

1.3 Quantum Solitons 
1.3.1 The particle theorists 

Solitons resemble extended particles, that is, they are non-dispersive localized 
packets of energy, even though they are solutions of non-linear wave equations. 
Elementary particles are localized packets of energy and are also believed to 
be solutions of some relativistic field theory. The particle theorists were thus 
highly motivated to find some quantum version of these classical solitons, that 
is, to quantize the solitons. 

It isn't immediately clear how to make the correspondence between a classical 
soliton and some extended particle state of a quantized theory, or between any 
classical field solution and its quantum analogue for that matter. To understand 
the difficulty, consider first the simple case of a point particle in a potential. 
Classically, this particle has some definite position and momentum with some 
particular path chosen by its initial conditions. Quantum mechanically, the 
picture changes entirely! No longer can we associate a particle with a definite 
position and momentum; instead, we must describe the particle probabilistically 
via a wavefunction 'ij}{x,t) giving the probability t)p to find the particle at 
point X and time t. How does one go from the soliton solution to some quantum 
wavefunction? 

Procedures for establishing this correspondence developed in the mid-70's were 
essentially a generalization of the semiclassical expansion of non-relativistic 
quantum mechanics. It was shown that not only could we associate a quantum 
soliton-particle with the classical solution, but also a series of excited states by 
quantizing fluctuations about the soliton 

For a soliton, we quantize its motion by defining conjugate position X and 
momentum P operators and imposing commutation relations. In the original 
field, however, there is an entire continuum of degrees of freedom that remain. 
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These are taken up by the quasiparticle excitations. 

The procedure is analogous to the quantization of a particle residing in a local 
minimum of the external potential (for example, a; = c in Figure n~B|) . This 
local minimum is not the global minimum, and hence is not the true ground 
state; however, there is a potential barrier blocking it from decaying to the true 
ground state. This is the same for a soliton excitation, or a vortex, which is 
higher in energy than the ground state, however, stable against decay. 

The quantization of the local minimum begins by assuming to zeroth order the 
classical solution, x = c. We expand the potential about this local minimum, 
finding quadratic behaviour to leading order, and proceed to quantize the per- 
turbativc excitations. Of course, a quadratic potential has simple harmonic 
excitations, so that the quantized solution can be envisioned as a hierarchy of 
simple harmonic excitations, centered, of course, about the classical minimum. 

For a soliton in field theory, the procedure is much the same. We begin by 
the classical solution, expanding the energy functional about it and quantiz- 
ing the leading order corrections. The simple harmonic analogous solutions are 
called mesons in quantum field theory, or quasiparticles in condensed matter. 
Of course, the mesons or quasiparticles also exist as excitations in the ground 
state, or vacuum state. Thus, quantization of the soliton is performed by ac- 
counting for the spectrum shift in the quasiparticle excitations and imposing 
commutation relations for the soliton position and momentum operators. 

For a good introduction on the quantization of solitons from the quantum field 
theorist's point of view, see the book of Rajaraman'^*' or the review articles of 
ColemaniS, or Rajaramanpi^. 

Recall, however, that the soliton is a spontaneously broken symmetry solution: 
in has chosen an arbitrary point in space about which to center. The Goldstone 
thcorc predicts a gapless boson mode restoring this broken symmetry. This 
causes divergences if we consider the next order semiclassical expansion of the 
quantized soliton, because of zero energy denominators that appear. 

An analogous situation for a simple particle is when the potential is completely 
flat. To all orders we find zero frequencies when expanding the potential. This 
is because all points are degenerate and the particle must randomly choose 
among them. In the quantum version, we find that the particle is no longer an 
eigenvalue of position at all, but rather of momentum, in the form of a plane 
wave. 

For the soliton, the Goldstone mode is dealt with in essentially the same way. 
For each broken symmetry, the quantized soliton has an associated momen- 
tum which is a good quantum number. For example, if the soliton exists in a 
translationally invariant system, we would find it has a well defined momentum 
in the quantized version. This, incidentally, provides a systematic method for 
calculating the mass of the soliton. 
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The general methods for separating the Goldstone mode involve introducing a 
collective coordinate for each broken symmetry>i2i22iiS. Since the original system 
doesn't depend on these coordinates, the final expanded energy functional can 
only depend on their conjugate momenta. 

The magnetic system of this thesis has a two dimensional translational symme- 
try broken by the introduction of a vortex. Thus, we promote the vortex center 
coordinates to collective coordinates to we obtain an effective action depending 

2 

only on the associated conjugate momentum via a particle-like ^ term. 



1.3.2 In condensed matter theory 

In condensed matter, we are more specifically interested in the physical con- 
sequences of the quantized solitons, as opposed to their mere existence and 
basic properties. Shortly after the quantum field theorists developed the soliton 
quantization methods, Krumhansl and Schrieffer showed that one dimen- 
sional quantized solitons could be treated exactly as elementary excitations, in 
addition to the ever-present quasiparticles. To explain, suppose we've quan- 
tized a soliton in a translationally invariant system (of length L with minimum 
length scale /). In the most general case, we would find, in addition to the 
regular Goldstone mode, a finite number of quasiparticle modes localized to the 
soliton, interpretable as soliton excited states, followed by the usual continuum 
of extended quasiparticle excitations. Krumhansl and Schrieffer show that the 
total internal energy of the system can be simplified to 

U=(^j- NbNl"'^ keT + Nl°' (e^ + hgT + {m ~ l)kBT^ (1.1) 

where Nb is the total number of localized quasiparticle states, including the 
translation symmetry-restoring Goldstone mode. 

This represents the internal energy of a system with ( j — NhNjf*'') quasiparticle 
modes and N^"* particles of rest energy S° each having ^ksT translational en- 
ergy and thermal energy /c^T for each of the Nt — l internal modes. The average 
number of particles Nl°^ forming the soliton is calculated using thermodynamic 
relations once we define a soliton chemical potential. See Currie et al.i for more 
details of the complete thermodynamic description of the soliton as an ideal gas. 

Quantum vortices were first considered by condensed matter theorists as early 
as the 1940's by Onsager**^ in superfluid helium. FeynmanJ^i developed further 
the idea of these vortex lines to explain the dissipation mechanisms for a rotating 
superfiuid and conjectured that they may also be responsible for the superfiuid 
to normal fluid phase transition. Unfortunately, in 3D the problem is essentially 
unsolved, so that no details of a vortex driven phase transition have yet been 
developed. 

In 2D, the problem is more tractable, and in the 1970's, Kosterlitz and Thou- 
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Figure 1.7: The two degenerate dimer states of trans-{Cll)N, polyacetylene. 

Ies9^ detailed a phase transition due to the proliferation of dislocations. The 
theory applies equally to vortices. Below the transition, the free energy is min- 
imized by maintaining the vortcx-antivortex pairs bound; however, raising the 
temperature to the transition, the gain in entropy by unbinding the pairs bal- 
ances the increase in energy. 

In two dimensions, the energy of a dislocation or vortex diverges logarithmically 
in the system surface area, 

E = Eo\n^ (1.2) 

where Aq ^ is the smallest area in the discrete system, where a denotes the 
lattice spacing. 

The entropy associated with the dislocation also depends logarithmically on the 
area since there are approximately A/Aq possible positions for it to center on, 

S^kgln^- (1-3) 

where /c^ is the Boltzmann constant. Since the energy and entropy depend 
on the size of the system in the same way, the free energy, F = E — TS, is 
dominated by the energy term at low temperatures so that the probability of an 
isolated dislocation in a large system is vanishingly small. At high temperatures, 
dislocations appear spontaneously as the entropy term takes over. The phase 
transition temperature can be roughly estimated as = E^/kB- 

In the late 1970's, very important new phenomena were discovered indepen- 
dently by the particle and condensed matter physicists. Jackiw and Rebbi-*' in 
considering the Dirac equation in the presence of a soliton found it had fermionic 
I states; while, Su, Schrieffer and Heeger— were studying kinks in a coupled 
electron-phonon model for the quasi- ID conducting polyacetylene and found a 
neutral spin i soliton state. 

Restricting ourselves to the polyacetylene system, consider a one dimensional 
system of electrons in a tight-binding model interacting linearly with the lattice 
coordinate displacements (essentially, coupling the electrons and phonons). The 
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Figure 1.8: The band structure of polyacetylene, gapped due to the electron- 
phonon interactions. Note the two isolated electron states in the 
gaps are only in the presence of a kink. 

Hamiltonian of this system is then 



where u„ and Pn are the lattice coordinate displacements and their conjugate 
momenta, characterized by mass m and stiffness constant K . The electrons are 
denoted by creation/annihilation operators c| ^ and Ci^s at site i with spin s, with 
hopping constant to and coupling constant a with the lattice displacements. 

The ground state of this system is doubly degenerate and spontaneously breaks 
reflection symmetry (this was predicted by Peierls*^ using mean-field approxi- 
mation for any non-zero electron-phonon coupling) . Figure 11.71 shows the two 
degenerate dimer states. As a consequence of the two-fold degeneracy, there ex- 
ist the kink and antikink topological solitons connecting the degenerate ground 
states (see Figure [131 — in actuality, the kink is spread over ~ 14a). 

Su et al.^ found that the kink had two states: a charged state, Q ~ ±e, with 
spin s = 0, and a neutral state with spin s = ±i. In addition, when the kink 
is in its neutral state, there is an s = electron state in the middle of the gap 
(see Figure [QH note there are two states, one localized to the kink, the other 
to the antikink) formed by pulling ^ a state per spin out of the Fermi sea. 

The polyacetylene study introduced to condensed matter physics what the par- 
ticle theorists independently introduced within a relativistic field theory: the 
existence of states with fractional charge. Although the i charge is obscured 
by the doubling of degrees of freedom due to spin, the zero energy state is 




N 



+ a 




(1.4) 
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Figure 1.9: A kink solution connecting the two degenerate dinier ground states, 
shown, left, on the linear polyacetylene chain, and, right, on the 
idealized chain with periodic boundary conditions. 

still formed by drawing half an electronic state (of each spin). Furthermore, 
the spin-charge relations are also unusual: charged solitons are spinless while 
neutral solitons carry spin ^. 

1.3.3 Superfluid He^ 

Returning our discussion specifically to vortices in condensed matter, quantum 
vortices were first proposed by Onsager*^ and developed more completely by 
Feynman— . A quantum vortex can be imagined as a regular fluid vortex with 
a cylindrical core shrunk down to atomic dimensions. The circulation of the 
vortex is quantized in units of h/m, where h is the Planck constant and m is 
the bare "^He mass. 

Describing the motion of superfluid vortices by making analogy to the motion of 
their parent fluid vortices was extremely successful. Early experiments by Hall 
and Vinen^^*^ found that if they applied an impulsive force setting a superfluid 
vortex into motion the vortex underwent helical motion (resembling that of an 
electron drifting in a magnetic field). In general, such a force arises always when 
a body with a flow circulation around it moves through a liquid or gas as in, for 
example. Figure IT. 101 

First noted in 1852 by Magnus when studying inaccuracies in the firing of cannon 
balls, the force responsible, named the Magnus force after its discoverer, can be 
explained in terms of the Bernoulli equation. The speed of the fluid is effectively 
lower on one side of the rotating body than the other (perpendicularly to the flow 
of the fluid, of course) so that the side with higher speed has lower pressure — 
thus the body experiences a force in that direction (see Figure ll.l()|l . The 
Magnus force in a superfluid is written 

Fm = PsK X (v - Vs) (1.5) 

where ps is the superfluid density, v is the vortex velocity and Vg is the asymp- 
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Figure 1.10: The equi-pressure lines of a fluid surrounding a rotating cylinder. 

The pressure differential top and bottom creates an upward force. 
The fluid flow is to the left. 

totic superfluid velocity (affected, of course, by the vortex presence). 

Hall and Vinen found the motion of their experimentally observed vortices could 
be explained with such a perpendicular Magnus force and an inertial mass of 
the order p^^, where p is the fluid density and ^ is the vortex radius. 

In addition, damping forces acting on the vortex were introduced with phe- 
nomenological parameters. The most general damping can act both longitudinal 
(as we are most accustomed to) and transverse to the vortex motion, expressible 
as 

= D(v„ - v) + X (v„ - v) (1.6) 

where v„ denotes the normal fluid velocity, whose exact definition might vary 
from one formalism to another. Note that the transverse damping term has 
the same behaviour of the Magnus force (with potentially an additional force 

OC V„ - Vs). 

Although this heuristic description is very successful in explaining observed 
phenomena, the microscopic derivation of the various parameters is far less 
successful. There is considerable disagreement, especially in calculations of the 
transverse dissipation parameter. 

An early calculation by lordanskii-^SiSi revealed a transverse damping force, later 
termed the lordanskii force, proportional to the normal fluid density 

Fj = PnK X (V - V„) (1.7) 

due to the scattering of phonons on the vortex. This entails an effective Magnus 
force with the superfluid density replaced by the total fluid density, pe p, plus 
additional forces proportional to v„ — v^. 

In the early 1990's, Thouless, Ao and NiuiSiSl (TAN) claimed that the transverse 
force was exactly the bare Magnus force of equation (|1.5|l . at all temperatures 
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while accounting for the scattering of phonons. The force on the vortex hne 
due to phonons is simply the variation of the phonon energy expectation with 
vortex position 



where fa denotes the occupation probability of the phonon state a. By ex- 
panding the phonon wavefunction to first order in vortex velocity using time- 
dependent perturbation theory, TAN were able to rewrite the force as an integral 
over the Berry phase associated with a closed loop around the vortex. Assum- 
ing no circulation in the normal fluid density, this reduces exactly to the zero 
temperature Magnus force. 

The transverse force on the vortex line can also be expressed as the commutator 
of the X and y components of the total momentum operator 



Applying Stokes' theorem, the integral over the cross-sectional area can be ex- 
pressed instead as a line integral about the boundary of the one particle density 
matrix. TAN argue that this boundary may be extended very far from the vor- 
tex core so that contributions from localized phonon states at the vortex core 
do not influence the transverse forced. 

In opposition to Thouless, Sonin— explained the transverse damping force via 
an analogous mechanism to the Aharonov-Bohn>i effect of an electron passing 
a double slit in the presence of a magnetic vector potential (though in regions 
of no magnetic field). The electrons passing in one slit relative to the other 
experience a phase shift due to the vector potential term, causing a horizontal 
shift in the observed interference pattern. However, this entails a momentum 
transfer from the magnetic field source, here a conducting coil, to the electrons, 
transverse to the double slit screen, and thus a transverse force acting on the 
coil. 

Similarly, quasiparticles passing above or below a moving vortex experience a 
relative Berry's phase shift*. A momentum transfer must occur between the 
vortex and quasiparticles, again, entailing a transverse damping force. 

Sonin calculated the effective transverse force exactly in the form 



so that the effective Magnus force is the regular Berry's phase result plus the 
lordanskii force. The normal fluid velocity here is in the vicinity of the vortex 
and may differ from the asymptotic velocity due to viscous dragging of the 
normal fluid by the vortex motion—. 

One apparent source of disagreement, first noted by Sonin, is that the vortex 
undergoes oscillatory motion due to the passage of phonon quasiparticles. The 




(1.8) 




(1.9) 



Ft {Ps + P7i)k X (v - v„) 



(1.10) 
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scattering calculations of Fetter— and Demircan et al^, which supported the 
TAN Berry's phase calculation, effectively held the vortex fixed by an external 
pinning potential, thereby nullifying the transverse damping force. 

The transverse dissipation is not the only source of controversy. The effective 
mass itself of the quantized vortex has not been agreed upon. Initial estimates 
are based on the inertial mass of the circulating fluid, essentially, pr^, with rg 
the radius of the vortex. In the quantum limit, the vortex radius shrinks down 
to atomic dimensions, or zero, so that the vortex mass tends to zero also. 

Alternatively, as suggested by Duan and Leggettii, the mass of the vortex must 
be proportional to 

cx ^ (1.11) 

^0 

where My is the vortex mass, is the stationary vortex energy, and vq is the 
velocity scale of the superfluid quasiparticles. This can be explained by purely 
dimensional arguments. 

For a quasi-2D vortex, however, the stationary vortex energy is log divergent in 
the system cross-sectional area, as in H1.2|l . suggesting the effective mass is also 
log divergent, much larger than the vanishing estimate made earlier. 

Clearly, the microscopic derivations of superfluid vortex dynamics has yet to 
firmly agreed upon. The variety of conflicting results suggests we re-examine 
the different methods used. Doing so in the simpler magnetic system is an aim of 
this thesis, though, unfortunately, a comprehensive study of the various methods 
could not entirely be undertaken. Rather, we calculate results here using regular 
perturbation theory, expanding in vortex velocity, and using Fcynman- Vernon 
influence functionals^. 

1.3.4 Magnetic vortices 

Magnetic systems have received much attention for their variety of applications 
and their lucrative potential-^, for example, in the market of magnetic memory. 
Vortices in magnetic systems are very easily observed and manipulated, for 
example using Brillouin light scatting— or magnetic force microscopy (MFM)^. 

Despite the ease of experimentally observing magnetic vortices, there have been 
relatively few microscopic derivations of the dynamics of vortices in magnetic 
systems. In fact, these derivations should be greatly simplified in a magnetic 
system; however, the resulting dynamics still possess many of the same strange 
aspects discussed with respect to superfluid vortices. 

A magnetic vortex experiences a force transverse to its velocity, the gyrotropic 
force. This force acts exactly in the same manner as the Magnus force, how- 
ever, has a different microscopic origin. It arises from a self induced Lorentz 
force, with the vortex vorticity acting as an analogous charge, while the out 
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of plane spins create an effective perpendicular magnetic field (this analogy is 
more fully developed in Chapter . Notably, this force is dependent on both 
topological indices (and is absent entirely for in-plane vortices for which the 
polarization is zero), as compared to the Magnus force dependent solely on the 
vortex circulation in a superfluid. 

There are interactions with quasiparticles that may alter the effective gyrotropic 
force. However, there have been no attempts to describe a transverse damping 
force in a magnetic system. In fact, all descriptions of dissipation in a vortex 
system have focussed on calculating an average energy dissipation rate or have 
been phenomenological (except for the work of Slonczewski^ which we describe 
in a moment). 

The earliest theoretical work on two dimensional magnetic systems with vortices 
are adaptations of the work of Thiel o^^i^^ . Thiele first introduced the gyrotropic 
force and dissipation dyadic acting on a magnetic domain wall in a three dimen- 
sional system. His dissipative force, however, was phenomenological employing 
a Gilbert damping parameter (the phenomenological damping parameter nor- 
mally introduced into the so-called Landau-Lifshitz equations governing the 
magnetization dynamics) . 

In the early 1980's, applying the work of Thiele, Huber^S^ and Nikiforov and 
Sonin^^ independently described the basic motion of a magnetic vortex. They 
calculated the gyrotropic force and phenomenological damping forces acting on 
a single vortex. 

Slonczewski— shortly thereafter considered perturbations about a moving vor- 
tex, deducing an effective mass tensor. A collection of vortices behave strongly 
coupled and the inertial energy is not diagonal but rather must be expressed as 
■^MijViVj where there is an implied double sum over the vortex indices i and j. 
He calculated the vortex dissipation via a frequency dependent imaginary mass 
term by studying the asymptotic behaviour of the lowest order vortex-magnon 
coupling. We will compare our dissipation results with those of Slonczewski in 
Chapter H 

Scattering phase shifts have been calculated for a variety of planar magnetic sys- 
tems^^''*^'^"'^. They were primarily interested in the thermodynamic behaviouii^ 
of such systems and searching for a vortex signature that could be measurable 
to verify a Kosterlitz-ThoulesSi24 transition. In fact, based on the modified spin 
correlations due to the presence of vortices, a central peak found in neutron- 
scattering experiments could be reproduced—. 

In a series of papers^^i^LSS, Mertens et. al. modeled numerically the motion of 
a vortex pair assuming various boundary conditions. The ensuing motion was 
best reproduced assuming an non-Newtonian equation of motion which included 
a third time derivative of the vortex position. 

We find just such a small third time derivative term in our influence functional 
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analysis. We compare our results with Mertens et. al. in section B. 3. 41 

However, this is a misapplication of the collective coordinate formalism: each 
collective coordinate is meant to replace a continuous symmetry broken by the 
vortex. In a planar system, a vortex breaks the two dimensional translational 
symmetry allowing the introduction of a two dimensional center coordinate only. 

There has been no work yet to find effective damping forces acting dynamically 
on a magnetic vortex. In this thesis, we calculate these forces assuming an 
averaged motion of the perturbing magnons. 



1.4 Easy plane insulating ferromagnet 

We study an insulating plane of spins, that is, fixed on their lattice sites, ferro- 
magnetically coupled, lying preferentially in the plane. The order parameter of 
the easy plane ferromagnet lies on the unit sphere but with an energy barrier 
at both the north and south poles. There arc hence topological solitons sponta- 
neously breaking the ground state symmetry, the continuous in-plane symmetry, 
and, at some energy cost to restore continuity, twisting out of plane to break 
the discrete up/down symmetry. There are also discontinuous vortices lying 
entirely in the plane as found in the XY model. 

We noted in the symmetry breaking discussion that a vortex lying in this or- 
der parameter space does not have topological stability. This however is for 
a completely degenerate sphere. Here, there is an energy barrier for paths to 
cross the two poles so that any homotopy of a vortex to a point would require 
passing a macroscopic number of spins through this energy barrier. The vor- 
tex thus has approximate topological stability, unless the anisotropy becomes 
vanishingly small. 

The energy of a general state {S^} of this lattice is 



where the indices extend over all lattice points in the 2D lattice. The first term 
is the exchange term and is approximated by including nearest neighbour inter- 
actions only, negative to ensure ferromagnetic coupling. 



where < i,j > denotes nearest neighbour pairs. For simplicity, we've assumed 
a constant exchange parameter J. The second term enforces the easy plane 
anisotropy, where K is the anisotropy parameter (for S > 1/2). 




(1.12) 
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Since we are interested in the low energy behaviour, we eliminate the short 
length scale fluctuations by describing the system in a continuum approxima- 
tion. Instead of a spin Si at site i, we now have a spin field S(r). Sums are 
replaced by integrals over space. For instance, the anisotropy term becomes 

and the exchange term becomes 

-1 Yl • S:'^ i E ^(S, - S,) . (S, - S,) ^ i /■ d\JiVSf 

<i,j> <i,j> 



where adding the constant 5^ terms doesn't affect the dynamics. Note that 
(VS)^ {VS^f + C^Syf + {VS^f. The redefined constants are given by 
J = J/2 and K = K/a^, noting that we use new dimensions for an anisotropy 
density. From here on, we drop the tildes and simply use J and K for the 
continuum versions of the exchange and anisotropy parameters. 

The Hamiltonian describing the system is then written 

H = s^J d\ {i{vef + sin^ e (^^(v^)^ - xj^ (i.is) 



where the spin field is expressed in angular coordinates, S = S* (sin cos 0, 
sin ^ sin 0, cos 6*). 

As explained in Appendix 1X1 and —ScosO are conjugate variables in the 
discrete lattice so that the Lagrangian can be expressed in the continuum limit, 
^ ^ J (Pr/a'^ where a is some lattice spacing length scale, 

^-^/ ^[-^os9.j,-^(^iVer+sin'e(^{V^r-^^^^ (1.14) 

where we've defined the speed scale c/r„ with c = SJa^ and the length scale 
= ^/J/2K. 



Using Hamilton's equations (|A.5|I or the Euler-Lagrange equation (jA.2p . we find 
the equations of motion 

c ot smO 
1 09 

-— =sin9y^e + 2cos9Ve -Vcj} (1.15) 
c at 



There are two families of elementary excitations: the pcrturbative spin waves, or 
magnons, and the vortices. The vortices have two forms: the so-called in-plane 
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solutions with polarization 0, and the out-of-plane solutions with polarization 
±1. The treatment in this thesis considers explicitly the out-of-plane solutions, 
however, setting the polarization to recovers the results for the in-plane solu- 
tions. The out of plane spin behaviour cannot be solved analytically; however, 
the core and far field asymptotic limits suffice for obtaining general results. 

The spin waves are small amplitude oscillations about the ferromagnetic ground 
state or about a vortex state, in both cases with an ungapped spectrum. The 
difference in the two spectra can be attributed to the vortex presence and yields 
an effective zero point energy to the quantized vortex. The equations of motion 
for the vacuum magnons are modified to the equations of motion of magnons in 
the presence of a vortex. The additional terms are interpreted as the magnon- 
vortex interaction terms. 

There is a one magnon coupling with the vortex velocity. Normally, considering 
a central system coupled to perturbative 'bath' modes, we find to lowest order 
a one magnon coupling with the vortex field. There is no such coupling here 
because the vortex is itself a minimum action solution of the same system in 
which the magnons arise. Thus, there are no first order variational terms. This 
assumes, however, that the vortex profile is unchanging in time. Allowing it 
to move about the system introduces a first order coupling between the vortex 
velocity and the magnons. 

There is also a two magnon coupling affecting the magnon energy with long 
range effects. This term scatters the magnon modes and hence alters their zero 
point energy. We attribute this shift instead to the quantized vortex state. This 
two magnon coupling has other dissipative effects and energy shifts that are not 
treated in this thesis. 

Wc first review the basic characteristics of the vacuum magnon modes and the 
vortex solutions. The gyrotropic and inter- vortex forces are found immediately 
by expanding the Lagrangian about a many vortex solution. 

We then examine the effects of the various couplings between magnons and vor- 
tices. The one magnon coupling can be interpreted as small vortex deformations 
when moving at velocity V or, alternatively, as a single magnon scattering event. 
The second order perturbation energy correction of this one magnon coupling 
goes as V'^ and is thus interpretable as an inertial energy, from which we can 
deduce an effective vortex mass. There is an additional imaginary energy shift, 
or a dissipation, from this coupling. 

The two magnon scattering term has a zero point energy shift and other magnon 
occupation dependent energy shifts. We do not retain higher order scattering 
terms, keeping only one magnon couplings, although they may indeed contribute 
more significantly to the vortex dissipation^i^*'^ . 

The dynamical effect of the one magnon coupling is examined fully in the 
Feynman- Vernon infiuence functional formalismii. The two sub-systems are as- 
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sumed initially non-interacting with the magnons in thermal equilibrium. They 
are thereafter allowed to interact, the magnons generally shifting out of equi- 
librium. The effect of the magnons is then averaged over by tracing out their 
degrees of freedom. This yields, in an averaged way, the effect of the magnons 
on the vortex motion. As found in perturbation theory, the one magnon cou- 
pling is responsible for two new terms in the vortex effective action: an inertial 
energy term and a damping force term. 

In addition to the usual longitudinal damping force, wc find a transverse damp- 
ing force reminiscent of the lordanskii force in superfluid helium. Such a term 
has not before been suggested in a magnetic system. The damping forces possess 
memory effects — that is, they depend on the previous motion of the vortices. 

For a collection of vortices, we find that their particle-like properties arc not 
independent. They have mixed inertial terms such as ^MijViVj and damping 
forces due to the motion of one vortex acting on another. 

Next, we review the basics of the two elementary excitations, first the magnons 
and after the vortices. 



Chapter 2 

Magnons 
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The plane of spins with easy plane anisotropy has a degenerate ground state. 
The spins are ferromagnetically coupled and thus prefer to align, however, they 
may choose to align along any direction in the plane — an example of sponta- 
neously broken symmetry. The Goldstone theorem predicts that there should 
then exist boson quasiparticle excitations that are not gapped. In this system, 
these Goldstone modes are the small amplitude, or perturbative, spin waves. 
When quantized, the excitations are termed magnons. 

The magnon spectrum in the easy plane ferromagnet is ungapped, however due 
to the hard axis, the spectrum is not simply the regular ferromagnet spin wave 
spectrum w{k) oc k^. Instead we find a spectrum with reduced density of states 
near w = 0. 

We begin by examining the small amplitude equations of motion satisfied by 
the magnons; thereby deriving the magnon spectrum and density of states. We 
calculate a few old results using spin path integrals as illustrative examples 
that we will need in later calculations. We derive the quantum propagator, a 
calculation following closely that of a simple harmonic oscillator. The quantum 
propagator is then manipulated to again reveal the magnon spectrum and, under 
a simple substitution, to yield the thermal equilibrium density matrix. 



2.1 Magnon equations of motion 

The magnons are the quasiparticle excitations of our system. As such, to de- 
scribe their motion and properties, we expand in small deviations about the 
ferromagnetic in-plane ground state 

^=e- 7r/2 

ip=cf> (2.1) 

where we've chosen the ground state </> = amongst the continuum of ground 
states without loss of generality. The complete system Lagrangian in terms of 
these perturbing variables ip and •& becomes 
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^™ = ^ / 5^ (^^ - i (-^VV - + I) ) (2.2) 

where J is the exchange constant and K is the anisotropy constant, a is the 
lattice spacing, c = SJa^ and 7-2 = 5^. The conjugate momentum is now 
Sd, the hnearized version of —Scos9. We essentiaUy expand the Lagrangian 
to second order perturbations to obtain a simple harmonic-like Lagrangian. 
Consequently, many calculations to come here mimic very closely those for a 
simple harmonic oscillator. 

Varying (|2.2|l with respect to (p and -d yields the magnon equations of motion 



Idif 

c dt 

1 d§ , 

=V' 

c dt 



(2.3) 



Alternatively, we could have linearized the system equations of motion, (|1.15|) . 
directly with identical results. 

The analysis proceeds in a plane-wave expansion. This system of equations can 
be solved by Fourier transforming so that V'^tp —k'^ip^. and V^i? —k^dk- 
Assuming harmonic time dependence, the eigenvalues of the equations of motion 
yield the magnon spectrum 

uj{k) = ckQ (2.4) 

where Q = ^k^ + The spectrum is not gapped (note the overall factor of 

fc), a reflection of the continuous degeneracy of ground states. However, the 
density of states goes as ^2^q2 remaining finite as a; ^ in comparison to the 
isotropic ferromagnet with density of states ^ diverging for zero frequency. The 
two systems are compared in Figure ITTl 

Alternatively, Fourier transforming the magnon Hamiltonian directly via 

m =J^ I d^fce-'^-^k (2.5) 
diagonalizes H„i = J ^S-dip — Cm to 

Hm - Y 1 ^ (fcVkV'-k + Q'^9k^?-k) (2.6) 
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Figure 2.1: A comparison of the easy plane magnon spectrum and density of 
states with the regular isotropic ferromagnet. 



To quantize the magnons, we impose the commutation relations between the 
conjugate variables (p]^ and Sdk 

[S^,, ^w] = -in{2ir r ^\ " > (2.7) 

We diagonalize the system now via the transformation to creation/anniilation 
operators 



normalized such that [ak,a|j.,] ~ (27r)^ ^ . Substituting for ip^ and in 

terms of ak and aj^. into the Fourier transformed Hamiltonian gives after some 
manipulation 



where ojk is again the magnon dispersion relation (|2.4(l . 

We interpret the operators and ak exactly as for the simple harmonic os- 
cillator creation/annihilation operators. The combination aj^ak is thus the 
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magnon number operator nk and the spectrum has energy ?iwk for each of 
the rifc magnons plus an additional zero-point energy ^huji^ for each wavevector 
k. 

Notice throughout that we associate factors of to the spacial and frequency 
integration measures to keep them dimensionless. This is consistent since the 
integrals replace sums appearing in the original discrete system. 



2.2 Quantum propagator 

The quantum propagator is an operator describing the time evolution of a quan- 
tum state. Although the vacuum propagator of the magnons is not needed for 
future calculations in this thesis, its calculation offers a simple application of 
spin path integration in our easy-plane ferromagnet. With only slight modifi- 
cations to this derivation, that is with the addition of a perturbing term, or 
forcing term, we obtain the quantum propagator for magnons in the presence 
of a vortex. We must save this calculation for later after we've derived the 
appropriate forcing term. 

Suppose initially we know the state of the system of magnons which can be 
represented in the tp basis. To find the state of the system at a later time, T, 

^i^, T)^ j difi'Kiip, T; 0)ij{cp\ 0) (2.10) 

where 

iHT 

K(^,r;^',0)= (^|exp —\^') (2.11) 

is the quantum propagator expressible as a path integral (see Appendix [C)l 

K{^,T-^',Q)^ j\[^{v,t),§{v,t)]e^p(^-j^ dtC„,[^,^]^ (2.12) 

and where Sm — Jq dtCm is the action with the Lagrangian Cm given in H2.2() . 

Before proceeding with the semiclassical approximation — here exact since we 
have no terms of higher order than quadratic — first Fourier transform to diag- 
onalize the problem in fc-space. Introducing the Fourier pairs of (p and i?, H2.5|l . 
the Lagrangian becomes 

Cm^sJ (^k^?-k - I (A:Vk^-k + g't?k^?-k)) (2.13) 

The path integration measure is now the product of these Fourier coefficients 
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I?b(r,t),i?(r,t)] ^ n^^kWd^kW 



Now to find the classical contribution, the equations of motion arising from this 
Lagrangian are 

iVSn=0 (2.14, 



■37 <:<3' ; V -"k 

The general solution with boundary conditions fk{0) — (p[^ and (pk{T) — ip]^ is 

t/^k ^ _ yk / sinwfc(r ~t) \ , f'k f sinwfe(r - 1) 
/ ^ sinwfcT V ^ cos ujkt J sinw^T V -^coswfe(T-i) 

(2.15) 

where Wk = cfcQ as before. 

Substituting the classical solution back into the action, after some simplification, 
yields the classical contribution to the action 

/a?(Pk Sk 
i2.r 2Qsin..r + ^•^^-) ^^^'^^■^ " '^-^'--^ ^'-''^ 

To evaluate the effect of quantum fluctuations, we solve the relevant Jacobi 
equation (adapted for a spin path integral as described in Appendix [2)) 

iVf'^O (.17) 



with initial conditions (/3(0) = and S'd(t) — 1. The determinant of the fluc- 
tuations is then given by ix{T) = iSQ/ksintOkT for each k. Combined with 
the prefactors in the path integration measure S/h, we find that the Gaussian 
integral over fiuctuations yields the prefactor 



/ ^ (2.18) 

]] 2TTihQ sin UkT ^ ' 

Assembling the various pieces, the propagator of the unperturbed magnons is 



Sk f a^cPk iSk 



(27r)2 2hQsmujkT 

(('^k'^'-k + </'k<y5-k) coaujkT ~ 2(^k</5-k) ^ (2-19) 

where ipk and (^k arc the Fourier components of the boundary functions ifir) 
and (p'{r). 
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2.2.1 Spectrum via tracing over the propagator 

By manipulating the propagator, we can recover the magnon spectrum. To 
explain, consider first the propagator of a single particle starting from position 
qq at time and going to position qx at time T 

iHT 

K{qT,T;qo,0)^{qT\exp —\qo) (2.20) 



Taking the trace of this operator, i.e. set qT = qo and integrate over the end- 
point Qo of the periodic orbit, we find 



iHT 

G{T) ^ I dqo{qo\exp ^ko) 



dqo ^ {qo\£,n) exp ^ 



h 



= ^exp^ 



iEnT 



(2.21) 



where {£,n} denote a complete orthonormal set of eigenstates of H. Using the 
normalization condition of these ^„ then yields the excitation spectrum of the 
Hamiltonian. 

The trace of the propagator H2.19(l thus provides another means to find the ex- 
citation spectrum. Set (pk = Vk ^^^'^ integrate over each Fourier coefficient 



// Sk ( iSk \ 
dVk\ hr-^7T'- ^^^P TFT- ^</?k</?-k (cos w^T - 1) 
^ y ImriQ smuJk-L \nQ smuJk-L I 

For ease of notation, assume that the product function applies to everything to 
the right of it, notably implying an integration over k within the exponential. 

Performing the Gaussian integrals 



r (T) = V\ / / ttTiQ sin uJkT 

' i-i-y 2nihQsmojkTy -iSk{cosujkT -1) 



and noting that coswfcT — 1 = — 2sin^ (^^^)' ^^^^ reduces to 
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1 



1 _ Q-iuJkT 



k 




k n=0 




(2.22) 



where {jik} denotes a set of integers n^. Thus, comparing with equation (|2.21|l . 
we find the excitation spectrum hujk{n + as expected. 

This method of recovering the excitation spectrum is of course only useful in 
the special case that G can be cast into this final form. Nonetheless, by taking 
the limit T — > 0, the ground state term dominates the summation so that we 
can always at least find the ground state energy. 

2.3 Thermal equilibrium density matrix 
2.3.1 Magnon density matrix 

A quantum state is represented by a wavefunction ■(/'(r, t). Generally, this state 
is a superposition of the system energy eigenstates, {^i}. For example 



The probability of finding the ith eigenstate upon measurement is and, by 
conservation of probability, = 1. This is a pure quantum state. Alterna- 

tively, a system may be a statistical mixture of eigenstates. In that case, the 
quantum state isn't expressible as in H2.23|l . but, rather, is described by a set of 
probabilities pi of finding the system in eigenstate upon measurement. 

We may have a pure quantum state describing the entire interacting system, 
which to some extent is the entire universe. Of course, we may then only be 
interested in a small subsystem within the whole. We wish to describe its 
quantum state only in terms of the subsystem coordinates. 

The density matrix is a notation for describing a quantum state, necessitated 
by statistical mixtures such as a thermal equilibrium state, or entangled states 
of two sub-systems for which each individual system must be described by a 
density matrix even though the complete system may be in a pure state. As the 




(2.23) 
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name implies, we express the quantum state by a matrix describing the density 
of the subsystem or mixture in terms of its eigenstates or coordinates. More 
specifically, for a pure state, the density matrix is 

Ptj = c^Cj (2.24) 

where the Ci are the coefficients in (|2.23f) . For a mixture, 

Ptj = ^ijPi (2.25) 

where pi are again the probabilities of finding the system is state . 

Supposing we have a pure quantum state, we can write the density matrix in 
the coordinate basis 

p(x,x')=^6(x)p.,c;(x') 

=^(x)V/(x') 

The vector x is broken into the coordinates of interest x and remaining coordi- 
nates q such that x — (x, q) . The reduced density matrix for the subsystem of 
interest is found by tracing out the uninteresting degrees of freedom 



/5(x, x') = J dqp(i, q; i', q) 



This effectively averages the effects of the external system. For a pure state, 
trp^ 1. It can be shown that trp^ is maximal when the ensemble is pure; for 
a mixed ensemble trp^ is a positive number less than one. 

A quantum system in thermal equilibrium has its eigenstates populated with 
probabilities given by the Boltzmann weighting factor e~^^^ . The thermal den- 
sity matrix in the coordinate basis is 

p(x,x') = 5]6(x)?*(x')e-''^' (2.26) 

i 

This should be normalized by the partition function ^ e^^^^; however, we will 
omit it for case of notation. 

But this form is extremely similar to the quantum propagator when also ex- 
pressed in this basis 

K{x, T- x', 0) = ^ e.(x)4* (xOe-*/"^^' 

i 

= (x|e^*/'"'^^'|x') 

Under the substitution T — > ^ihf3, in fact, we recover the thermal density 
matrix, though unnormalized. See Appendix lA. II for formal details of this sub- 
stitution. 
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Noting this imaginary time correspondence between the quantum propaga- 
tor and the density matrix in thermal equihbrium, we make the substitution 
T — > —ihP in the quantum propagator 



Sk / f a'^(fk iSk 



(27r)2 2QsmwkT 
{{'fiW-k + Vk^-k) COS WfeT - 2(/Jk</j'_k) ^ 
to obtain the thermal equilibrium density matrix 



, / Sk / f a'^cPk Sk 

Pi'P, <P ) - liy 2nnQsmhniOkl3 \ J (27r)2 2gsinh?kjfe/3 

{{fW-k + </5kV'-k) coshhtOkfi - 2v3kV-k) ^ (2.27) 

This corresponds to the magnons being excited such that a state with energy E/- 
is measured with probability weighting given by the Boltzmann factor, e~*'^^'= . 



2.4 Summary 

In summary, the easy-plane magnons perturbing the vacuum ground state have 
the spectrum 

u!{k) = ckQ 



where Q = ^k"^ + ^, c = SJa^ and r1 = 



We calculated the real time propagator of these magnons and consequently, 
making use of the imaginary time path integral of the density matrix, also the 
thermal equilibrium density matrix. The propagator is extremely similar to that 
of a simple harmonic oscillator. In fact, under the substitution ^ — > mw the 
magnon propagator becomes identical to that of the simple harmonic oscillator. 



Next, we examine the vortex excitations. 
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Vortices 



The easy plane fcrromagnet admits two families of elementary excitations. In 
the last chapter, we reviewed the perturbative excitations, the magnons. Now, 
we review the other elementary excitations, the non-perturbative vortices. 

Although the out-of-plane spin behaviour cannot be described analytically, we 
present the asymptotic behaviour which is sufficient for getting leading order 
results. By superposing many vortex solutions, we expand the action to reveal 
an inter-vortex Coulomb-like force. The analogy is complete with the corre- 
spondence of 47reogi with electronic charge in Coulombs. 

The dynamic term "pf?" in the action is rc-cxprcsscd describing a gyrotropic 
force (analogous to the Lorentz force) or, alternatively, as an effective dynamic 
term in terms of vortex coordinates, P • X, where the momentum term is a 
vector potential. This is analogous to a charge in a magnetic field for which the 
momentum is modified by the magnetic field vector potential. In this formalism, 
the corresponding vector potential describes an effective perpendicular magnetic 
field B = £^PiZ. 

We briefly present different possible two- vortex motions: depending on the rela- 
tive sign ofpiQi, the pair execute parallel motion (for opposite signs) or co-orbital 
motion (for like signs). This basic motion is perturbed by introducing an iner- 
tial mass term. Finally, the zero point energy shift of the two niagnon coupling 
is examined in a Born approximation. This approximation is found to be suf- 
ficient for the continuum of magnons; however, there exist translation modes 
localized to the vortex core. These will be reconsidered in the next chapter 
using collective coordinates. 

The system has two symmetries: a continuous in-plane symmetry and a discrete 
up-down symmetry. The vortices are thus characterized by two topological 
indices, the vorticity q = ±1, ±2, . . ., sometimes also called the winding number, 
and the polarization p = 0, ±1. 

The p = vortices are often separately considered, termed the in-plane vortices, 
while the p ^ solutions are called the out-of-plane vortices. This separation, 
however, is unnecessary: allowing p = 0,±1 in the following treatment recovers 
the proper results for both types of solutions. 
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Being non-perturbative solutions, the vortices satisfy the full, non-linear, equa- 
tions of motion of the easy plane ferromagnet. Derived from the system La- 
grangian 

^ = ^ / 5^ (-<^ ^ - i ((V^)^ + sin^ ((V0)2 - 1) ) ) (3.1) 
the equations of motion are 



i^=-^+cos.(V0)^-^eos^ 
1 36 

-— ^sm9V^e + 2coseve -Vcj) (3.2) 
c ot 

where J is the exchange constant and K is the anisotropy constant, a is the 
lattice spacing, c = SJa^ and = 

The in-plane vortex can be described analytically. The spin configuration of this 
solution has (ji^ = q£, + S and 9^ — 0. The parameter q is called the vorticity of 
the vortex, and (5 is a phase that has little importance on the vortex dynamics^ . 
We can solve for its energy within our continuum approximation, requiring both 
an infrared and ultraviolet cutoff, 

E = S^ [ Srii^c^^f = 5Wln— (3.3) 
12 a 



where Rs is the radial size of the system and a is a lower cutoff, the lattice spac- 
ing, required since the system is actually discrete (making r — > unphysical). 
Note that this energy is independent of where the vortex center is within the 
circular integration region. 

The out-of-plane solution is also characterized by its polarization; that is, the 
direction (up/down) that the spins twist out-of-plane. The spin configuration 
has the same polar angle dependence, 4>v = qS, + <5, while the out-of-plane spin 
angle cannot be solved for analytically. The asymptotic behaviour is 




where ci and C2 are free constants that can be set by imposing appropriate 
continuity conditions. 

Figure ini shows the spin configuration of two simple out-of-plane vortices. This 
solution has the same leading order energy as the in-plane solution 

^Note that this broken continuous symmetry entails the existence of gapless boson modes: 
the magnons. 
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<p<P cOi=o<i=i<t] cOi=o 



Figure 3.1: Vortex spin configuration: left, a vortex with q = —1; right, a vortex 
with q — I. 



Ey ^ S^jT:q^\n— (3.5) 
a 

Core corrections to the energy are finite and hence negligible in comparison to 
this log divergent contribution. In fact, in most that follows, the core will be 
ignored since it usually offers a finite contribution next to a log divergent one. 
A notable exception is the gyrotropic force that depends on the core behaviour 
via the core polarization. This is a differentiating feature of magnetic vortex 
dynamics from that of classical fluid or superfluid vortices where the analogous 
Magnus force depends only on the vortex circulation, the fluid analogue to the 
magnetic vorticity. 

The motion of the in-plane vortex undergoes many of the same corrections. In 
fact, with the substitution p the treatment here reduces to that of an in- 
plane vortex. The gyrotropic force disappears, however, all other forces and 
correction are polarization independent. 



3.1 Force between vortices 

Consider two vortices of vorticity qi and q2 and polarization pi and p2 well 
separated so that the only distortion in their profiles can be assumed to lie in 
the region between the two where their profiles are entirely in the plane. The 
spins in this middle ground are aligned in the plane with angle determined 
by the sum of spin angles (see Figure rO|l given by each vortex independently 

012 = qix(Xi) + q2x{X2) (3.6) 



The out-of-plane component of the spin can be neglected here since we've as- 
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Figure 3.2: Two vortex spin configurations. Left, two vortices with q—\\ right, 
vortices with q~\ and q — —\\ both with no relative phase shift. 



sumed that the vortex cores are widely separated and each core gives only a 
small correction. 

The energy of the two vortex system is 

Ei2^^l d\ [{ye,,f + sin^ {Vcj^.^f + (3.7) 

which, except for regions within radius of each vortex core, is dominated by 
the (V(/)i2)^ term. Thus, neglecting core terms, the energy becomes 




As an illuminating trick to evaluating this integral, note that 

( '^^^^ I 92^\ _ ( glXl 92X2 \ 

But £ — + ^^f^ is just the electric field generated by a pair of point 

charges, Airegqi at Xi and 47reo(72 at X2, in two-dimensional electrostatics using 
SI units. The electrostatic energy, including the divergent self energies, of this 
configuration is exactly 

W=^h,^ + jLi^E + Alr,^ (3.9) 
27reo Tv 27reo Ty Treo ry 

where X12 is the vector from vortex 1 to vortex 2. Alternatively—, we can 
express the electrostatic energy as the integral oi ^£^. Thus, upon comparison, 
the energy of the two vortex system is 
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Figure 3.3: Intervortex forces: top, two vortices of opposite vorticity attract; 
bottom, two vortices with same sense vorticity repel. 



E,2 = S^Jn {q'(ln^+ ql In^ + 2q,q2 In ^ ) (3.10) 



Similarly, for a collection of n vortices, with cores widely separated, the spin 
field pattern is 



tot =^gix(Xi) 

i=l 
n 

= ^^,(r-Xi)«0 (3.11) 



^tot 

i=l 



Following the same analogy to electrostatics as before, we find the energy of the 
collection of vortices is now 

Etot =S^JwJ2 <ll In — + 2S^Jt, ^ ^g. In ^ (3.12) 
The force acting on vortex j due to vortex i, separated by distance 
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Fij = — V Xij Etot 

(3.13) 

Ay 

where Xy is a unit vector pointing from the center of vortex i to the center 
of vortex j. Thus, if the two vortices have the same sense, or the same sign 
vorticities qi and q^, the force is repulsive, and conversely, for opposite senses 
the force is attractive. Note, since in this approximation there is no interaction 
between the two vortex cores, the direction of the spins out of the plane at the 
cores — the polarization — is irrelevant. 



3.2 The gyrotropic force and the vortex 
momentum 

3.2.1 The gyrotropic force 

The vortex is a stationary solution of the system. If we assume that it now 
moves at a small velocity X, for the moment with no deformation to the vortex 
profile, the pq action term, called the Berry's phase in a spin system, is no longer 
vanishing. The Berry's phase, lob, is a phase accumulated by the changing spin 
field 



u}B= dt —ScosOct) (3.14) 



a 



Considering a single spin, we can interpret this phase geometrically as the solid 
angle swept out by the motion mapped onto the spin sphere. This is clear when 
we make the change of variable 



lub = S 



J d(j3Cose = S J doj' 



where diu' is the area increment on the unit sphere. Refer to Figure 

We treat this term as a potential and calculate the corresponding force acting 
on the vortex. Let the vortex profile move as a function of r — X(t), where X(t) 
is the center coordinate of the vortex. The Berry's phase term in the Lagrangian 
becomes 

f d^r ■ f d^r ■ 

—S / cos9y — S / — ^X • V(/)t, cos6'„ 

7 J 

The gyrotropic force arising from this term is found by varying it with respect 
to the center coordinate of the vortexiSiS^, without the usual negative sign since 
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Figure 3.4: The spin path mapped onto the unit sphere. The area traced out 
by its motion gives the Berry's phase. 



we take the term from the Lagrangian, 



f (fr ■ 

Fgyro =S'9x / ■ ^'Pv COS9v 

= -S J ^ V (x • cos ej^ (3.15) 

But the integrand is strictly a function of r — X so that 9x —V, where V is 
understood to be with respect to r. Note that V'^t/*,, = and thus 

V ■ Vcpy cos ej^ = ■ V cos 

Using the cross-product relation A x (B x C) = (A • C) B — (A • B) C, we find 

(^X- V(/)„) Vcos^„ = S {Wcos0v X V0^)xX-(^X- Vcos6'^) 

where now both terms on the right are integrable. Consider the first term, not- 
ing that 

(Vcos6' X V0 ) _ d cos 9yd(j>y d cos 9 ydcj)^ _ d {cos 0y,4>y) 
^ ^ ^ dx dy dy dx d{x,y) 

Clearly, since V cos 9y and V(pv both lie entirely in the plane, the z component 
is the only non-zero component. The first integral becomes 

/cPr S f 2tt Spq 
cos9v X Vdy = — dcos9ydd)y z = — z 
J 
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A Fgyro A 

^ A ^ f ^ 



gyro 



Figure 3.5: The gyrotropic force: left, a vortex with p — 1 and q = —1 travehng 
to the right experiences an upward force; right, a vortex with p = 1 
and (7=1 travehng to the right experiences a downward force. Note 
z is defined out of the page. 



where p is the polarization of the vortex core and q is the vorticity of the vortex. 

For the second integral, consider axes and xj_ parallel and perpendicular to 

X, where the second is aligned such that z x X = xj^. In polar coordinates 
defined for this frame, the integral can be written 

~~2 J (i^ ■ y cos 9y^ \7 (fly — J drdxX cos x(— sin x, cos x) 

where we decompose x = (~sinx,cosx) into the {xi\,x±) basis. Evaluating 
this gives 

" J / • Vcos6'„) V(/)„ = TTpqXxjL = X X (3.16) 

The gyrotropic force is then 

Fgyro = X X (3.17) 

Note, this result differs by a factor of 2 from that of Hubert using the formalism 
of Thiele, found for a magnetic domain wall^^. Thiele's starting point for the 
kinetic term was ^ 

— cos 9y + — cos 6y (3.18) 

at 

which is exactly twice our starting point that includes only the first of these two 
terms. 

Notice that the gyrotropic force is derivable from the equivalent Lagrangian 
term 
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L-gyro = ^ X X • Z - X X Z • X (3.19) 



3.2.2 The vortex momentum 

The gyrotropic force can be written in the suggestive form 

dP 

^gyro = -J^ (3.20) 

where Pgyro is a momentum term from the equivalent Lagrangian term H3.19|l 
written in the form P • X 



P...o = -^Xxz (3.21) 



We now examine a direct evaluation of the vortex momentum as given in a 
general field theory by the operator— 

P = - y (frn{r, t)V^{r, t) (3.22) 

where tt is the conjugate momentum density to the field variable (p (the tilde's 
are there to differentiate the field variable here to the azimuthal angle (j) used 
previously). This operator is chosen because it is the infinitesimal generator of 
spatial translations, eg. 



0(r + 5v) = (j)[v) + 64>(r) = 0(r) + V0(r) • Sr 
5^{r) = {6r ■ P, 0(r)} = V(^(r) • 6r 

where we use the Poisson bracket here as defined in equation IA.8I (note there 
we used q for and p for tt). 

For the magnetic vortex, this gives the momentum expression 

S 

d^r^ cos 6l.„V0„ (3.23) 

Before attempting to evaluate this expression, first note that the 1/r behaviour 
in V(/)^ is balanced by the r in the integration measure so that the integrand is 
nowhere divergent. 



If we blindly set the vortex at the origin of the integration region, the Xr di- 
rection of the integrand sums to by symmetry, there being no other angular 
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dependence. The integral is non-zero, however, if we displace the vortex by X 
from the origin. 

To evaluate this integral note that 

V0„ = -qz X Vln|r-X| (3.24) 

Considering the momentum integral one component at a time, first the y com- 
ponent 



/ 



A.™M.ln|r-X|=-/rf.d!,4™«.,(ln,. 



rX 



2 



drdxrdr COS 9^ cos Xrf ■ X 
irqpX 



where X is the x component of X. We expanded the In above and truncated the 
series to 0{l/r). This is in keeping with the r — * behaviour noted in the orig- 
inal integral. Of course, for r oo the integrand decays to zero exponentially 
as before. 

After the analogous treatment for y, we find the momentum is exactly the Fgyro 
describing the gyrotropic force 

Poi/ro = o ^ ^ (3.25) 

What does it mean exactly to have a momentum that is speed independent 
and coordinate dependent? Isn't this extremely bizarre? Recalling the problem 
of a charged particle in a magnetic field, the momentum of such a particle is 
modified by the presence of the magnetic field according to^ 

p ^ p - -A (3.26) 

c 



where A is the vector potential describing the magnetic field B = V x A, e is 
the electric charge and c here is the speed of light. 

For the magnetic vortex, this momentum term must also correspond to a vector 
potential term. Completing the analogy, using Aireoq as charge as in section 
13.11 replacing the speed of light by the speed of magnons SJa'^/ry, we find an 
effective perpendicular magnetic field B = -^^pz. 

To further explore this interpretation, we expect the gyro-momentum to be 
gauge dependent. That is, we should be able to rewrite the vector potential 



A ^ A + V,/(r) 



(3.27) 
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for any continuous function /(r), changing the momentum expression Pgj/roj 
however, and stiU describe the same physical system. 

Considering this gauge change in reverse, we use the gauge freedom of the 
Berry's phase. The Berry's phase is written in a general basis 

^B = j dt(frA{n)n (3.28) 

where A is a unit magnetic monopole vector potential. We change the gauge of 
this vector potential A via 

A^A + V^/ (3.29) 



where / is a general function of n. 

The momentum of the magnetic vortex is altered by noting the correspondence 



^ j,gyro — J" d vA^^j.^^^ 



The Berry's phase gauge change shifts the momentum definition according to 

/S 
<fr— cos9yVr(t)v + Vj^./VrOi (3.30) 

But 

Vfj./Vrf^^ = Vr/(fl) = -Vx/(fl) 



since tl = r2(r — X). Thus, the additional term to the vortex momentum 
becomes 

- / d^rVc^JvA = ^^j d'rf = Vxi^(X) (3.31) 

where FpC) = J d^rf is now some general function of X. 

A continuous function -F'(X) can always be expressed as the integral over another 
function /(r,X). Thus, the gauge freedom in the Berry's phase allows exactly 
the necessary gauge freedom in the gyrotropic momentum term, supporting our 
vector potential interpretation. 
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Figure 3.6: Sequenced photographs of a pair of fluid vortices with same sense 
vorticity. Photos were taken at 2 second intervalsi^. 



3.3 Motion of vortex pairs 

Consider the motion of a pair of vortices, separated enough that the cores do 
not significantly interact, with polarization pi and vorticity qi, i — 1,2. The 
motion so far is dictated by the balance of the inter- vortex and gyrotropic forces 
acting on each vortex 



Tvpiqiz X Xi =0 

7rp2g2Z X ±2 =0 (3.32) 



or taking the cross-product of each equation with z 

Y 2S'^Jpiq2 , , 
Xi = X2 = -2 z X (X2 - Xi) 

P2qi ^12 

In the case ^192 = Xi = — X2 and the vortices move on a common circular 
orbit, keeping separation X12 with angular frequency w = 4S Jpig2 ^jjg^-g ^ > Q 
denotes counter-clockwise rotation. 

For the opposite case, Piq2 ~ ^P2qi, we have Xi = X2 and the vortex pair move 
with a common velocity on parallel lines (upward for P2qi > and downward 
for ^291 < 0). 

In this approximation, the dynamics of the vortex pair is identical to the anal- 
ogous motion of a pair of fluid vortices. Referring to Figures H3.6|l and l|3.7|l . 
a pair of fluid vortices with the same direction circulation move in a common 
circular motion while a pair of opposite circulation move along parallel paths. 
There is the notable difference, of course, that here the type of motion is dictated 
by the products pq rather than just q as in regular fluid dynamics. 



2TrS Jqiq2 



^12 



2TTS^Jqiq2 



^12 



(Xi ~ X2) 



(X2 - Xi) 
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Figure 3.7: Sequenced photographs of a pair of fluid vortices with opposite sense 
vorticity. Photos were taken at 4 minute intervals^. 



3.4 Vortex mass 

Up to now, we've assumed that the vortex profile is rigid when in motion. In 
fact, the profile is modified linearly in X. 

Assuming that the vortex moves at small velocity X, expand about the rigid 
vortex profile 

=0l°)(r-X) + 0W 
9, (r-X) + 0« 

Substituting this into the equations of motion, (|3.2|) . to first order in X, making 
use of the zeroth order equations of motion, these reduce to 




(3.33) 



+ sin 20(0) V(/.,« -V^o) 
--± ■ V0(") = sin0(o)v20(i) + 2cose' (v0[°) • V^l^) + VO^^^ ■ V(/.(°)) 

Using the asymptotic expansion of the stationary vortex (|3.4|) , the asymptotic 
forms of the profile perturbation become, keeping only the dominant terms in 
the r — > and r — > oo limits 



0(1) = 



-^^^cosx, r->0 
^cosx, r — > OO 



"^"^ sinx, r —> Q 

Xc2rl'^ e-/-- _ _ _ (3-34) 



where Ci and C2 are free parameters in the unperturbed asymptotic form H3.4|l . 
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Substituting these asymptotic expressions into the energy integral, we find en- 
ergy terms that are quadratic to lowest order in X (the linear terms integrate 
to zero by symmetry) interpretable as a ^MyX"^ kinetic term: 

E^E(^^+Eill + E(^^ (3.35) 



where Ecore accounts for the r = O-.r^ and E^^ accounts for the remaining 
r — r„..oo. Evaluated, 



/o 2 2 1 2 2 5 



p(l) ^ / o 2 2 _ ^ 2 ' 



2 2 _ ^ 4 2 



8Ja4 V ' " 6 ^ 12 



2 / 15 ^2 „2^2 



Assuming an energy correction of the form 



AE = ^AUX^ (3.37) 
the leading term describing the vortex mass is deduced as 

M. = ^ln^. (3.38) 

Note that this mass is, in fact, identical to the mass estimate suggested by Duan 
and Leggettii based on purely dimensional arguments, M„ = j^^^ry- 



3.5 Quantization of magnetic vortices 

Quantum fluctuations in a system introduce a zero-point energy. In the previous 
chapter, we quantized the magnons finding this zero-point energy to be ^ \h,ujk, 
summed over the entire fc-spectrum. In the presence of a magnetic vortex, the 
magnon spectrum is shifted. Since we prefer to have a consistent definition of 
the magnons and vortices, the shift in the zero-point energy of the magnons is 
associated with the quantized vortex. 

Quantization of a magnetic vortex involves quantizing the small variations about 
it and examining how the energy of these modes shift from the analogous modes 
in the absence of a vortex—. See Appendix IB . II for more details. 



Expanding 6 and <j) about a vortex, 9 = Oy+'d and = + (/s, in the non-linear 
equations of motion H3.2|l . yields the linearized equations in d and 
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sin6l,„VV + 2cos6'„ {VO., • + V?9 • V0„) (3.39) 

These are very similar to the vacuum magnon equations of motion with the 
addition of a few perturbing terms. Notably, these additional terms all decay 
away the vortex core and will be treated in a Born approximation (applicability 
of this approximation is discussed at the end of the next section). 

Alternativelyi22i22i^j we could expand as 



smOv dip 

c 'at 

1 dd 
c~dt 



=9, + d 

M J. , ^ 



smt 



(3.40) 



yielding the linearized equations in d and Lp 



-^=~ V^-d - cos26'„ (4 - (V0„)^ 1 + 2 cos e.,Vip ■ V0„ 
c ot \ry 

=VV + cos2 6l^, - (y(j)yf ] ip+iVe^f ip + 2 cos 0yVi9-V4>y 



or, equivalently, in the more symmetric form 



C ot ^ ' X'^ Ok, 



where a; = radial derivatives are now with respect to x, and 

Vx(x) = ((V0„)^ - l) cos26l^, 

Y2{x) - ((V0,)' - l) cos' 0„ - (V^„)' (3.42) 

This form is particularly suitable for examining the core effects and searching 
for possible bound modes. 
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3.5.1 Phase shifts in the Born approximation 

The perturbing terms are locahzed to the vortex and can be treated in a Born 
approximation. The dominant scattering term decays as whereas the re- 
maining terms, neglected in the foUowing treatment, die exponentially. The 
error introduced by neglecting these terms will be discussed in the final analysis 
at the end of this section. The magnon equations of motion are modified to 



c dt V 
c~di 



~ =VV (3.43) 



The perturbation treatment is most straightforward in a single variable. Elimi- 
nating the ?? variable, we have 

^^'^ -VV+f^-^)v^^ (3.44) 



Note the additional term ^ modifying the vacuum equations of motion of the 
magnons (|2.3() . The Born approximation is applied using the standard partial- 
wave analysis from scattering theory'''^. 

Consider the orthonormal basis functions such that V^^k —k^^M and 
assume harmonic time dependence. We expand Lp in this basis 

(p = ^Ck'e^'^'=V (3.45) 

k' 

where to zeroth order we've assumed 

" [0, otherwise. 

The zeroth order terms simply reduce to the vacuum equations of motion. The 
first order terms are 

- E - - E fc'^Q'^^VCk^ + ^^k (3.47) 

k'^k ' k'^k 

where we've cancelled the common e"^*"* factor. Recall = fc^ + p-. Multi- 
plying by i^^„ and integrating over space, enforcing orthonormality of the {^k}, 
we find an expression for the first order coefficients 
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Substituting for the unperturbed magnon spectrum (|2.4|) and using plane waves 
for the orthonormal basis, the first order correction to (p is 

(27r)2 fc2Q2 _ ^,2Q/2 J ^,2^' 

First, integrating over the polar angle 0k' from to tt, we obtain 

(1) 1 dk' k' Hi'\k\v-V'\)+H^^\k\v-V'\) /■ ,2 , -.k'.r'?!^! .U.r' 

The ±ie are chosen to displace the poles so as to pick the outgoing wave (the 
plus is for the -ffp^^ integral, the minus for the h'^^ integral). Considering the 

f 2) 

i?Q integral, there are poles in the complex k' plane at k' = ±(fc + ie'). Noting 
the asymptotic behaviour 



H[^\kr) ^ J^e-^('='-('+^)^) (3.49) 
V irkr 

we close the contour about the positive imaginary axis for the Hq integral, to 
pick up the k' ^ k + ie pole with residue 2{k'^+Q'i) ^o^^ ~ '''D- ^^'^ integral 
over Hq^^ is just the complex conjugate (c.c) of that over Hq'^^ and hence follows 
immediately. 

Thus we have 

^^^^ W - ^^^^^ I d\'H^,'\k\r - r'D^e--' + c.c. 
We now expand the Hankcl functions according to the identity 

OQ 

H^''^\k\r - r'l) = ^ Ji{kr')H[^'^\kr)e'^^^-^'^ (3.50) 

/— oo 

if r > r', an allowable assumption if we only want the wavefunction correction 
for asymptotic r, and (f) {(j>') is the polar angle of r (r'). Expand the plane wave 
as 



' = i"'Jm{kr')e'""f'' (3.51) 



After integration over (f)' (giving a factor 27r(5;m), the wavefunction correction is 

2 ^ r 2 

^'"'^'^ = Aik^ + Q2) T. ^'Hrikr)e''J dr'±jjfikr')+c.c. (3.52) 
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Recall we assumed the unperturbed solution was a plane wave, expandable 
according to (|3.51() . so that the entire solution can be written, up to first order, 

^(r) = Y: [Mkr) + ( 4(J%2) gr^(fcO j dr'^-jfikr') + c.c.J 

Comparing this with the sum of an incoming and outgoing wave 

i [e'^'H\^\kr) + e-'''^'i/f )(fcr)) ^ Ji{kr) - H['^\kr) + H'l^\kr) 

^ -e-'^' cos (kr + -)- - A/ ) as r -> oo (3.53) 



Trfcr V ' 2'2 



gives for the phase shift of the 1*'* order wave 

^' (3.54) 

These phase shifts perturb the magnon wavevector k — /c«'="** — A/, and, hence, 
the magnon spectrum LOk- For proper counting of the total energy shift, first 
discretize k by fixing the boundary conditions of the wavefunction at r = Rs so 
that 



Tin = fc„i?5 = fcr"-Rs - AK^n) (3.55) 



Notice that asymptotically we have a cosine wavefunction as opposed to a plane 
wave as described by Rajaraman— . Letting the system size tend to infinity then 

The zero point energy shift, given by the change in the zero point energy of 
the small oscillation modes when the vortex is present as compared to those in 
vacuum, is then 

Ai; i ^ nSujk (3.56) 



2 

k,l 



But k = fc"'^'^" - A 
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duj{k) Ai{k) 



dk Rs 

so that 



A ^ ^ f dk du} ^ . ,„ ^. 



Substituting for Ai{k) from equation H3.54|l . noting that ^ = 



4 7 fc2 + g2 Q 



/ — — OO 



4 7 Q r^, 

oo 

where we've used that ^ Jf{kr') — 1. Note that the radial integral is cut off 

/ — — oo 

by the vortex core size. This is because the perturbing term changes behaviour 
drastically in the core so that our analysis cannot be extended there. The k 
integral can be evaluated noting that 

The result is ultraviolet divergent so that we must impose a cutoff of 1/a, 
physically reasonable if we recall that a is the lattice spacing of the discrete 
lattice. Finally, the energy shift in the presence of a magnetic vortex is 

A^; . In ^ f - ^ In :i±^M±^\ (3.58) 

4 y 2a?r^ 2rl " / 

This zero-point energy shift, due to the presence of the vortex, is associated to 
the quantized vortex rather than the magnons^*'. Thus, Ai? is the zero-point 
energy of the vortex. Note that the interaction actually decreases the quantum 
energy of the vortex-magnon system. 

We can examine the error in neglecting the exponentially decaying terms by 
replacing the r^/r' behaviour by exp(— r/r^). Essentially, this would replace 
the log divergence in the final result with unity. Hence, in comparison with the 
main contribution, these exponentially decaying terms are negligible. 
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The Born approximation amounts to the substitution of 

(</)/|[/|V'^) ^ (</>/|t/|</'.) (3.59) 

where and '0 denote the unperturbed and modified waveforms, respectively. 
In general, the validity of the Born approximation depends on how much the 
waveforms differ in the region of the scattering potential**^. In our case, the Born 
approximation indicates that the two wavefunctions in fact differ to first order 
by equation (|3.52l) which is proportional to the predicted phase shifts. This is 
circular reasoning; however, in the case of those quasiparticle modes delocalized 
over the system, we expect the waveform not to change significantly. On the 
other hand, there are quasiparticles that become trapped by the vortex center. 
Clearly, for these modes, the wavefunctions are drastically modified in the vortex 
core, where the scattering potential is greatest, so that a Born approximation is 
invalid. We examine these bound modes in the next section to show how they 
result from the translational symmetry broken by the vortex solution. 



3.5.2 Bound modes 

As pointed out by Ivanov et. al^^iS^t^, the short range interactions neglected 
in H3.43|l can drastically alter the behaviour of certain modes. The symmetric 
perturbing equations, H3.41|l . are more suitable for exploring the core region. 
Assume a solution of the form 

I? =f{x) cos(mx + ojt + ip) 

If =g{x) sin(mx + ujt + -ip) (3.60) 

Substituting this into 13.41|l yields equations for / and g 

/ ^ 1 d vn? V { ) \ f ( + 

\ dx-^ X dx x"^ J \ c J 

+ 2-"^2(a;) U= + 2 / (3.61) 

\dx x dx x / \ c X / 

recalling that x — and that Vi (x) and V2 (x) are defined in (|3.42|) . 
For a; = 0, there exist exact solutions for m = 0, ±1 



q sin 9v 



(3.62) 
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For |m| > 1, the asymptotic behaviour of the modes is entirely unbounded so 
that the vortex center has not greatly shifted the magnon wavef unctions and 
the Born approximation applied in the previous section should be valid. 

Consider first the m = result. Combining the unperturbed vortex profile with 
this result (recall the normalization of the perturbations as in H3.40|l 'l 

9 =9^ (3.63) 

where 5x is the coefficient of the linearized solution. We find that it corresponds 
simply to the freedom of uniform rotation in the a;y-plane. 

Similarly, consider the m = ±1 solutions 

(j) =qx — - — sin{mx — V') 
r 

9 =9y + m6r9'y cos(mx — ijj) 

with Sr as the coefficient of the linearized m = ±1 solution. But note that the 
additional contributions can be re-expressed as 

=qx + V0„ • mSr 

9 =9y + V9^ ■ m5v (3.64) 

where 5v is now a vector of magnitude 5r in the direction defined by the polar 
angle tp (see Figure . Thus, these two modes represent infinitesimal motion 
along ±(5r (the sign chosen by the sign of to). 

Clearly, these bound modes are inadequately treated using the Born approxima- 
tion and must be treated separately somehow. Ivanov et. al.^2i224^ attempted 
to calculate the phase shifts of these modes separately and to subsequently use 
them to describe the angular and translational motion of the magnetic vortices. 
Alternatively, however, one can treat the problem using collective coordinates 
(see Appendix IB. 1.21 for more details) conveniently separating these so-called 
zero modes and treating the remaining modes in a Born approximation. 

In the next chapter, we expand the interactions of the magnetic vortex with 
the environment magnons using collective coordinates. Using a path integral 
formalism, we separate the degrees of freedom of the vortex motion from those of 
the environment and proceed to integrate out these modes yielding the effective 
dynamics of the vortex. 
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Figure 3.8: The directions relevant to a small translation of the vortex along i5r. 
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Vortex dynamics 



We now have all the background to interact the vortices and magnons. Using 
a variety of techniques, we examine the effects of couplings between the two 
systems to the vortex energy and dynamics. In the previous chapter, we already 
saw how a modification in the magnon spectrum can be interpreted as a quantum 
energy shift associated with the vortex. 

First, using regular perturbation theory, we examine the one magnon coupling 
with the vortex velocity giving rise to an inertial mass and a dissipation rate 
of a moving vortex. We also examine the long range two magnon coupling in 
this language, finding almost immediately the zero point energy shift that in 
the previous chapter required calculating all magnon phase shifts. 

The effective vortex dynamics are derived by finding the time evolution of the 
vortex-magnon density matrix and tracing over the magnon modes. We use 
the Feynman- Vernon formalism, describing the density matrix with path inte- 
grals. We again deduce the vortex inertial mass, in agreement with perturba- 
tion results. The vortex motion is again dissipative; however, we find the vor- 
tex damping forces explicitly and characterize the associated fluctuating forces. 
Generalization to a collection of vortices is carried out. 

In addition to the previously derived gyrotropic and inter-vortex forces, we 
derive microscopically vortex damping forces. We introduce for the first time in 
such a magnetic system a transverse damping force, analogous to the lordanskii 
force acting on a vortex in a superfluid. These are derived from the action terms 
found in the vortex density matrix propagator (|4.91|) and have corresponding 
fluctuating forces with correlations given by H4.92|l . 

Alternatively, we consider decomposing the motion in a Bessel function basis, 
{Jrn{kX{t))e'^"^'^^}, to obtain Brownian motion for the components with an 
effective action given by (|4.96|l and corresponding fluctuating force correlations 

Km . 
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4.1 Vortex-magnon interaction terms 



We work with the complete non-hnear Lagrangian for our magnetic system 
(frf- c ( , , . T „ 1 



I _ -0COS0-- (V^)2+sin2 (V<^)2-_) ) ) (4.1) 



Expanding the Lagrangian density about the vortex profile via 9 — 0^ + •& and 
(f> = (f>jj + ip we find the following terms in the integrand 



— + (^) (- cos 0^ + sin 0^d)- — (( V6'„)^ + 2W0y ■ Vi? + 
(sin^ 0^ + sm29y + cos 26*^ d^) (^(W(f>^)^ " ^ "^^^"^ ' ^ 

The zeroth order terms in Lp and 'd simply give the vortex action; the first order 
terms give § multiplied by the equation of motion and (/? multiplied by 
the equation of motion and thus are zero, except, notably, the one magnon 
dynamic term 

g . 

^(j)ysm9yi} (4.2) 
Finally, the remaining two magnon terms are 



a 



S 5'2j , _^ _^ _^ _j, 

— V3sin6'„ I? — UV'df +8111^ 0y{V Iff + 2sin26'„V0„ • Vip 1} 

J4 \ 



+ COS20, ((V0,„)2-l)i?2] (4 3) 



Minimizing these action terms, we find the perturbed equations of motion sim- 
ilar to 



sin0„ d^p , sine, ^ _ v^^? - cos20„ i\ ~ (V0,)^ ) i) + sin2e„V(^ • V<^„ 



c dt c \r. 
1 d-d 

-— :=sin6l.„V^</7 + 2cos6'„ (V6'„ • V(^ + Vt? • (4.4) 
c at 

Define the vortex profile relative to the center coordinate^ X 



t There is no need to add a collective coordinate reflecting the rotational symmetry of the 
problem since this is actually just a restatement of the 2-dimensional translational freedom, 
already entirely taken care of in the 2-dimensional center coordinate. 
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(4.5) 



The center coordinates play the role of the collective coordinates in this system, 
introduced to account for the continuous translational symmetry broken by the 
vortex. They are elevated to operators. 

Focussing on the one magnon perturbative term, (|4.2|) . expanding in terms of 
the collective coordinates, we find 



In the previous chapter, this term perturbed the vortex profile under motion, 
introducing a finite vortex mass. 



4.2 Perturbation theory results 
4.2.1 Vortex mass revisited 

Consider the one magnon coupling (|4.6I) . This term can be considered as 
a perturbing term of the vortex profile under motion or, alternatively, as a 
vortex- magnon coupling. Fourier transforming according to H2.5|l . now with 
I? = 'i?(r — X), we can rewrite the coupling as 



where the r integration has been shifted to move the vortex center coordinates 
into the exponential. Expanding Xvr — Xvk + Xkr and noting that 



-^<pv sin -d 



(4.6) 



a- 






the coupling term becomes 
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S f <Pka? /dVe''^'^™''^'='-sin^ sin x«fc cos Xfer 



-.2 



27r)2 " r 

(27r)2 ^'''"'"^^k y drX sinXvkJiikr) sinOy 



2TTiSq f (Pka'^ -ik x 
2'KiSq f (Pko? _j,^.x <j X • Xfc 



k- 



(27r)2 *^ fcQr„ 

where we approximate sin 0„ w 1 — e^''^''" , which has the right asymptotic 
behaviour for r — s- and r — > oo. Expressing i9k in terms of creation and 
annihilation operators given in 12.8(1 . we finaUy obtain 



X^^^t _ (4.7) 



(27r)2 a2 y 2kQ Qr^ 

or expressed a httle more symmetrically 



(27r)2 a2 Y 2fcQ Qr^, 



— zk-X t zk-X \ ^ A o\ 

e cik - e a_k) (4.8) 



Note the similarity of the coupling here to that of the polaron problem discussed, 
for example, by Feynman-'^^. 

In first order perturbation theory, this coupling provides no energy shift since it 
necessarily changes the number of magnons between the initial and final state. 

In second order perturbation theory, we consider the diagram shown in Figure 
14.11 corresponding to the emission and re-absorption of a virtual magnon. The 
energy shift provided by this diagram, which with foresight we call E„iass, 



(X • Xk) 



(27r)2 a^kQ^rl V V (P - Kkf /2Af„ + hckQ - P2 /2M„ 
However, so far the vortex has no inertial energy, P'^ /2M 0, and 

d^ka^ 1 2Tr^hSq^ / • n 2 



E 



(27r)2 hckQ a'^kQ^rl 



(X • Xk) (4.9) 



The integration over the polar angle of k contributes a factor tt. We expand the 
k and Q dependence in partial fractions 

1 rl krl k 



kQ^rl k Q2 g4 
The radial integral is evaluated as 
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P-k 



Figure 4.1: Lowest order contributing diagram for the first order vortex-magnon 
coupling term. 



2Ja4 J k'^Q^rl 
2Ja4 j \k Q2 Q4 

where we've imposed both an upper and lower cutoff, with a the lattice spacing 
and Rs the system size. 

Thus, identifying this as a ^M^X'^ inertial term, we find a vortex mass of 

My = In , ^ - " 7^ , o 4-11 

in agreement to leading order with the analysis of section 13.41 The r„ 
-\/a^ + replacement corrects the r„ — )■ limiting behaviour. 



Mass tensor of a collection of vortices 

We can easily generalize this result to a collection of vortices in this formalism. 
Recall that the n- vortex superposed solution is given by 



1=1 
n 

^tot (4.12) 

1=1 

so that the one magnon coupling becomes 
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Figure 4.2: Definition of angles for evaluation of off-diagonal mass terms. 
The second order energy correction is now 

= E / "'"^jM^, ^ ^ (X, ^ (4,14) 

The diagonal terms evaluate exactly as above. The off-diagonal terms can be 
evaluated noting that 

J dxknj sinxfet,! sinxfct,2 exp{ikrij cosxfcnj 

= \ J '^Xkr,, (cos(xtoi ~Xkv2) +cos(xfei,i -l-Xfei>2))exp(iA:r„ cosxfcr„ ) 

where = jX^ — If we assume the various angles are defined as in Figure 

^4.21 then Xkvi \kv2 \v\V2 \kv\ ~^ \kv2 \v\rij ~^ \v2rij ~^ ^Xfcrij • 

Noting by symmetry that the smxknj terms integrate to zero, we have 

j dxkr.j sinxfc^i sinxfc«2 exp{ikrij cosxknj) 

^\ j '^^kn, (cosx,;it,2 +cos(xt,ir,, + Xv2r,,) coa2xkr,,) exp{iknj COS Xkn^) 
= TrcosXviV2Jo{krij) -f tt cos(xt,ir,,- + Xv2ri,)J2{krij) 
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Next we perform the integrals over k, noting that 



dk =ln- 
dk- 



1/B.s k 



_J2{krij) _ 1 







k 2 
and rewrite 

COs{Xvirij + Xv2r,j) = (Xi • eA)(Xj • Ga) - (X,j X Ga) ' (Xj X Ga) 

Finally the energy correction term becomes 

n 

Emass = 2^^^^'^^^ (4-15) 

where Mij is the n-vortex mass tensor given by 



"tI^I -(X.xe,,).(X, xe,,)j, t ^ j; (4.16) 



where = jxhix^- This result is in agreement with that of Slonczewski— . 
Slonczewski's calculation follows very closely that of section 



4.2.2 Radiation of magnons 

In the previous section, we calculated the vortex inertial energy using second 
order perturbation theory. However, we only used the principle part of the 
integral. When evaluating the integral giving the second order perturbative 
energy shiftii^, to be careful in the divergent region Ef — > Ei, symbolically, we 
should write 

AE, = y^Mlfl_ (4.17) 

and then take the limit e ^ 0. But 

1 X ie 

X + ie + x^ + 
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The imaginary part approaches a (5-function as e ^ since 



J — ( 



So then 



principle value ( — ) — i'iTS{x) 



x + te \x 

An imaginary part to the energy shift creates a decaying exponential factor in 
the time-dependent wave function 



and is hence interpreted as dissipation. The factor of 2 is there so that the 
probability |V'P decays as e 

The rate of decay due to magnon emission is thus given symbolically by 
If we assume an initial state of no magnons the rate equation becomes 

Again, initially we have no inertial term, simplifying 7 to 

. = ^-/|^^(x.&)'M"C) (4.21) 

Rewriting the (5-function as 

and the rate of dissipation becomes 



222 



hJa^ 



(4.22) 
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Slonczewski— calculated microscopically a dissipation rate by extending the 
simple results of section to include retardation effects of spin waves. How- 
ever, to evaluate the far region perturbations, he assumes small frequencies 
which give a log divergent frequency dependent dissipation. 

In our analysis, we used the same vortex-magnon coupling, although, accounting 
only for the k = contribution. We find the full spectrum contributions in 
section [4.4. II and will return to this comparison then. 



Second order radiative corrections 

Interestingly, the rate of emission using the finite mass of the vortex calculated 
in section l4.2.1l is actually considerably more important. This is reasonable since 
when we assume a finite mass, the vortex can lose kinetic energy by emitting a 
finite energy magnon. 

Thus, with a finite My this time. 



2tt 



Note that the vector potential component of the vortex momentum won't be 
included because in perturbation theory we assume that the vortex momentum 
is changed by changing speed, not position. 

Let be a solution of the delta function condition as a function of emission 
angle, Xk- The delta function can then be rewritten 



(, (P cM„\V f, fP cM„^^^^ 



k- — cosxk T- - 77T-p fcx - ^ cosxk 



(4.24) 



where we've approximated Q ~ l/r^, a reasonable approximation assuming 
small vortex velocities. Changing variables within the delta function to express 
it as 5{k — k^), the integral becomes 



(X-Xk) 



Snq^r, /•^2^ A MJjk - k^) 



fcx-(fcosxk-fe^) 



(4.25) 



Substituting for fc^, 
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Figure 4.3: The dissipation rate from perturbation theory; first assuming infi- 
nite mass and then adding corrections due to finite mass. 



My6{k - k^) 



f cosxk 



' hPa^ 



sm Xk 



cos Xk- fr 



(4.26) 



\ RscMy P2 _ 1 



where P = To evaluate this last integral, an infrared cut-off had to be 

imposed: k„un = 

The discontinuity at P = 1 occurs when the vortex attains the minimum energy 
to overcome the "semi-gap" formed by the Q — y^?^ + l/rf factor in the energy 
spectrum. 

Note that this dissipation is in addition to that calculated in the previous section. 
We didn't get both contributions here because we left out the fc = solution of 
the (5-function H4.24|l . See Figure H31 for a plot of these two contributions. 
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4.2.3 Zero point energy 

Consider next the two magnon couplings, (|4.3|l . arising from expanding the 
Hamiltonian about a stable vortex. Separate out the terms corresponding to 
the magnon Lagrangian expanded about a vacuum solution 



and interpret those remaining as an interaction Hamiltonian 



Hint J ^y~ip{l- siney)!} + |( - cos^ 9^ {Vipf + 2sin26'i, V(^^ • V<^ -d 
+ cos20„ (V0,)2^2_2£^^2U (4 27) 



The only long range interaction term in Hint (i-e. that doesn't decay expo- 
nentially) is the sin^ 9y portion of the fourth term. Fourier transforming the "d 
factors according to (|2.5|l , this long range term becomes 



= -ScaqnJdr J J 

(4.28) 



This integral over r diverges in the short range. However, the original term 
actually changes sign as r — > so that the analysis is invalid into the core region 
anyway and must be cut-off. We define a form factor J-{k) = e"*'''"^/ dr 
where k = k — k' . 

Expressing this term in the language of quantized magnons (see section I2.1|l 
gives a Fourier transformed version 
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TiSJq^a^TT f (Pk cPk I kn 



r d^k d^K 

J (27r)2 (271)2-^^'^ 



fiSJq^a^-n f d^k d'^n ^, ^ I kn 



f d^k d^K 



(4.29) 



The a and terms above correspond, respectively, to the case of two magnons 
being created in opposite directions, two magnons incoming from opposite di- 
rections being annihilated, a magnon given a momentum boost of k, and the 
last combination removes momentum k from an existing magnon. The last term 
gives the zero point energy shift 



AE = - 



hcq^n R 
In — 



d'^k k 




2rl 



In 



(4.30) 



as found before in section HT^l 



4.3 Vortex influence functional 



In this section, we develop the effective dynamics of the magnetic vortex using 
path integration. The temperature is introduced by assuming as an initial con- 
dition that the magnons are in thermal equilibrium. They are of course allowed 
to evolve out of equilibrium when interactions with the vortex are considered. 

Populating the magnons at a temperature r, we have a density matrix describ- 
ing them given by equation H2.27|l . To describe the effective dynamics of the 
reduced density matrix for the vortex, we trace out the magnon degrees of free- 
dom from the full vortex-magnon density matrix using the Feynman- Vernon 
influence functional formalism Before proceeding, consider the simple case of 
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a central coordinate x{t) coupled to a bath of simple harmonic oscillators ri(t) 
with frequencies Wj. This introduces the influence functional formalism and the 
interpretation of results for our own magnetic system. 

Separate the Lagrangian describing a coordinate x{t) coupled linearly to a set 
of harmonic oscillators as 

C = C4x{t)] + Crin] + Cint[x{t),nit)] (4.31) 

where £j.[a:;(t)] describes subsystem x(t), Cr[ri\ describes the environmental 
modes and Cini[x{t),ri{t)] describes the couplings between the two systems. 

Assume a general Lagrangian £3;[a;(<)] for the central coordinate, a simple har- 
monic Lagrangian in r, 

^r[n] = Y.rf + Y''^ (4-32) 

i 

and for the interacting Lagrangian, assume linear couplings 

jCinMt), nit)] = ^ Cix{t)ri{t) (4.33) 



Generally, the dynamics of the two subsystems become entangled which is con- 
veniently described within the density matrix formalism. 

The density matrix of the complete system in operator form evolves from initial 
state p(0) according to 

iHT iHT 
p(T) = exp —p{0) exp (4.34) 

Alternatively, in the coordinate representation, 

p{x, n, T; y, Qi, 0) ={x, ri\p{T)\y, qi) 

= J dx'dy'dr'^dq'^{x,ri\exp—-^—\x,r^ (4.35) 

iHT 

X {x', r'M^W, exp ^|t/, qi) 

Expanding each propagator as a path integral, noting 

(a^,r,|exp —\x',r'i)= j V[x{t)] j P[r,(0] exp -5[a;(i), r,(t)] 

iHT 1'^ i 

(a;,ri|exp —\x\r',)^ j ^W)] J ^ V[qi{t)]exp--S[y{t),qi{t)] 
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the density matrix at time T becomes 



p(x,n-y,q,;T)= j dx' dy' dr[dq[ I V[x{t)] j V[n{t)]e^p^-S[x{t),r,{t)] 



n 



y 



y Qi 



x{x',r'MO)\y',<ll) I T^W)]l 2?feW]exp--5[y(i),ft(t)] 



However, suppose we're only interested in the dynamics of the subsystem 
regardless of the specific behaviour of the harmonic oscillator subsystems. To 
eliminate these variables, perform the trace over the {ri\ variables to obtain the 
so-called reduced density operator 



p{x;y;T)^ / dr, / dx' dy' dr'^^ / V[x{t)] / V[n{t)]ei^^^S[x{t),n{t)] 



y 

y' -'I'i 



X {x',r'MO)\y',q^ I V[y{t)] I V[q,{t)]e^p --S[y{t), q,{t)] 



Assume that the t = density matrix is separable in the two subsystems, i.e. 
that they are initially disentangled and 

p{x,ri]y,q,]Q) = p^{x,y]Q)pr{ri,qi]Q) (4.36) 



Further, assume that the simple harmonic oscillators are initially in thermal 
equilibrium so that Pr{ri, qi;t = 0) is given by 



TT / muji muj / 2 2\ 

Prin, q.;Q) = 11 ^ exp - ^^--^ {{r, + ) cosh hu^.P - 2nq. 

The reduced density matrix is then expressible as 

y;t^T) = J dx' J dy'J{x, y, T; x' , y' , Q)p^{x' , y' , 0) (4.37) 

where 



J(a;,y,T;x',y',0) = / V[x(t)] f V[y{t)]ex^^Us,[x{t)]-SMmW).v{i)] 

(4.38) 

is the propagator for the density operator and 
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F[x{t),y{t)]^ J dTMdq[pM,q[,0) V[n{t)] j^^ V[q,{t)] (4.39) 

X exp J {Sr[n{t)] + S^nt[n{t),x{t)] - Sr[q^(t)] - 5,„t (i), 2/(i)]) 

is the influence functionali^i. Evaluating this for the central coordinate x{t) 
coupled linearly to a set of environmental modes described by simple harmonic 
oscillators with spectrum oji (t) 



1 /"^ /"* 

F[x,y]—exp~— / dt ds {x{t) — y{t)) {a{t — s)x{s) — a* {t — s)y{s)) 
" Jo JQ 

(4.40) 

where 

■^-^ / / N 2 cos iOiit — s)\ 
ait -s)^J2 <^^P - + . a 1 (4-41) 

where Ci are the linear coupling parameters. 
4.3.1 Quantum Brownian motion 

Caldeira and Leggett^ interpret the influence functional result as quantized 
damped dynamics. The problem of quantizing Brownian motion was not entirely 
understood. Their idea of coupling to a bath of oscillators to achieve Brownian 
motion (which, of course, from there is easily quantizable) was one of many 
proposed in the 1980's and 90's. The classical equation of motion for Brownian 
motion, the Langevin equation, is 

mx + rix + V'{x) = F{t) (4.42) 



where m is the mass of the particle, 77 is a damping constant, V{x) is the po- 
tential acting on the particle and F{t) is the fluctuating force. This force obeys 



(Fit)) =0 

{F{t)F{t')) ^2r]kTS{t - t') (4.43) 



where ( ) denote statistical averaging. 
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With such a force, the propagator of the density matrix of system x is given by 

J{x,y,t;x',y\0)^ J V[x\V[y\exp'- [s[x\ - S[y] + j ^ dT{x{T) - y{T))F{T)^ 

Assuming that the fluctuating force F{t) has the probabiUty distribution func- 
tional P[F{t)\ of difi^erent histories F{t), the averaged density matrix propaga- 
tor becomes 

J{x, 2/, t- x\ y\ 0) = y V{x\V{y\V{F\ P{F{t)\ exp ^ (S{x\ - S{y\ (4.44) 

+ y*dr(ar(r)-t/(T))F(T)) 

We can perform the path integration over F(t) if we assume P[F(t)] is a Gaus- 
sian distribution, yielding 

J(x, y, t; x', y', 0) ^ J V[x\D[y]exp'- {S[x] - S[y]) (4.45) 

X exp - ^ / / dTds{x{T) - y{T))A{T - s){x{s) - y{s)) 
n Jo Jo 

where A{t — s) is the correlation of forces, (F(r)F(s)). 
The real exponentiated term in the influence functional is 



exp-i j j dTds{x{T) - y{T))aR{T - s){x{s) - y{s)) 



(4.46) 



where 

- ^) = E ^ ^ - (4-47) 

i 

where Cj denotes the coupling coefficient to the ith environmental mode. As- 
suming instead a continuum of k states with density 

^-(-)^» = {o,^ ' ltd ^'-''^ 

the influence functional result becomes in a high temperature limit (coth 

2kT \ 

hanir - s) = {F{t)F{s)) = 2r?fcT ''"^^^ (4.49) 

7r(r — s) 
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which tends to H4.43|l m the hmit il oo. 

The imaginary phase term in the influence functional is manipulated to give 
an frequency shift which renormalizes the external potential. In addition to 
this, there is a new action term corresponding to a damping force 



Note that the forward and backward paths are interacting so that the new 
effective action is coupled in x{t) and y{t). 

The relaxation constant 7 is 

= — (4.51) 

where the damping constant rj is dependent on the density of states of the 
environmental modes. For our treatment where the environmental modes are 
magnons, we know explicitly the magnon density of states, going as ^.-i'^qi 
(recall Lo{k) = ckQ) rather than lo'^ as assumed above, so that our analysis does 
not simplify to a frequency independent damping function. 

Castro Neto and Caldeira™ consider the problem of a central coordinate coupled 
linearly to a set of oscillators; however, as opposed to the Caldeira and Leggett^ 
problem, the central system, X(t), is a solution in the same medium as the set of 
oscillators. Hence, as in our problem, there is no linear coupling with position, 
but instead, we find a linear coupling CiriX{t) between oscillators {n} and the 



They simplify their results by assuming a Born approximation. Although they 
lose the resulting frequency dependent motion, they do find that the damping 
coefficients and correlation integrals now possess memory effects. We will discuss 
these issues after results have been simplified for our vortex-magnon system. 

4.3.2 Semiclassical solution of perturbed magnons 

Before evaluating the influence functional, we flrst need the propagator of the 
magnon system perturbed by the vortex presence. The effect of important per- 
turbing terms have been discussed already using perturbation theory. The one 
magnon coupling endows the vortex with an effective mass and makes the vor- 
tex motion dissipative by radiating magnons. The leading two magnon coupling 
provides an overall zero-point energy shift to the vortex-magnon system that is 
associated to the quantized vortex. Although the two magnon couplings, or in- 
deed any of the many magnon couplings, may give more signiflcant dissipation, 
we neglect these contributions in this treatment. In the influence functional, the 
forward and backward paths have cancelling zero point energy shifts and hence 
we will ignore entirely the many magnon couplings. 




(4.50) 



velocity. 
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Treat the disturbance of a magnetic vortex centered at X(t) with vorticity q 
and polarization p by the magnons via the one magnon couphng, (|4.6() . 

We must evaluate the propagator for the system of magnons, again in the ip 
basis 

iHt 1"^ i f'^ 

((^|exp — —\ip')= J V[ip]V[d]e^p-J^ dtiC^+Cnt) (4.53) 

where X(t) is considered now an externally controlled parameter. 
Introduce the plane wave decomposition l|2.5|l so that the action becomes 



~ j ^e-^'^-'X- V</>.sin0,i?k) (4.54) 

The equations of motion are modified by a force term, that, for simplicity, we 
denote as 



/k[X] = y ^e^-X- V./)„sin0, 



= ^|;^e*- (4.55) 



and become 



-i cQ' )\K )~ \ -./k[X] 



(4.56) 



The solution with boundary conditions </3k(0) = ^p'^ and V3k(r) = (pk is 



'/'k ^ c!scoswfe(T- s)/k[X] / sinwfet 

i?k / ~ sinwfcT V ^coswfei 

ft 



Vk + /o cosa;fcs/k[X] / sincjfe(r - t) 



sincjfeT V cos uJk{T -t) 



(4.57) 
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Substituting this solution into the action gives the classical contribution 

Sk 

S'rL+^nt{k) ^^Q^^—f {{^W-k + COS^fcT - 2<^k¥'-k (4-58) 

-2(^k/ dtcosujktf^k[X{t)] + 2ip'i, [ dtcost^fe(T-t)/_k[X(i)] 
Jo Jo 

+2^ dtj^ dscosLjfc(T-t)coscofcs/k[X(t)]/_k[X(s)] j 



The quantum fluctuations introduce a pre-factor given by solving the relevant 
Jacobi equation 



m ^Q^ J\ ^(t) 



(4.59) 



with initial conditions (pk{Q) = and S'd{0) = 1. The determinant is given as 
i(p(T). Combining the pre-factor (det)"'^ with path integration measure factors, 
give the overall result, (j2.18|l . 



Sk 



2TTihQ sin ujkT 



(4.60) 



The final propagator is 



i{Hm+int)t. n TT ^k if o^d'^k i 



2TrihQ sinuJkT h J (27r)2 

(4.61) 



4.3.3 Evaluating the influence functional 

Substituting the semiclassical solutions to the two path integrals and for the 
thermal equilibrium density matrix, the problem is reduced to three regular 
gaussian integrals, ignoring pre-factors, which cancel anyway after all integrals 
when the density matrix is properly normalized. 
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1. 



(4.62) 



exp 
exp 



2hQ sinh hcukP 
iSk 



2hQ sin ujkT 



( ('/'k'/'-k + '/'k'/'-k) cos ujkT - 2ipkip'_y. 

2(^k / dtcosuoktf-v[Mt)] + 2^'^ ( dicosc^fe(T-<)/_k[X(t)] 
Jo Jo 

2/ rft / dscosc^fe(T-i)cosc^fes/k[X(i)]/_k[X(s)]) 
Jo Jo ^- 



exp 



2?iQ smuJkT 

cT 



( (<^k'^-k + </'k<<5-k) coswkT - 2(^k'/'-k 



2(/5k / dtcosujktf-i,[Y{t)]+2<f'^^ f dtcoswk(T-t)/_k[Y(i)] 
Jo Jo 

2/ dt f dscosujk{T~-t)cosujksMY{t)]f-u[Y{s)]) 
Jo Jo ' . 



where the J^j^ appUes to everything (and hence imphes integrals over k within 
exponentials). 

Performing these integrals mimics very closely the calculations for the analogous 
problem of a central coordinate x(€) coupled linearly to the position coordinate 
of a system of simple harmonic oscillatorsi^. In fact, the final expression is the 
same, with the same substitution niLo — > Sk/Q found earlier in evaluating the 
magnon propagator (see section [T^ . 



J^[X,Y] =exp-i j dtj^ ds(/k[X(0] -/k[Y(t)]) (4.63) 

(a,(t - s)/_k[X(s)] - alit - s)/-k[Y(5)]) 



where 



_ ,) _ ^ i^-^^.^t-.) 2cos^fc(^-.) 



(4.64) 



The propagator of the density operator can be written as 
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J(X,Y,T;X',Y',0) = / 2?X / PXexp^ 

r dt £ (/k[X(0] - /k[Y(i)]) - s) (/-k[X(s)] + /-k[Y(s)]) 
- i (/k[X(i)] - /k[Y(t)]) a^it - s) (/-k[X(s)] - /-k[Y(.)]) 



5^[X]-5^[Y]- (4.65) 



(27r)2 



where ct^{t — s) and a^(f — s) are the real and imaginary parts of ak{t — s) 

n , , Sk . , , hcOkP 

a^{t - s) =— cos(Jk{t - s) coth 

Sk 

aiit- s) =- ^sinujkit- s) (4.66) 

Ordinarily, we would extract a spectral function J{uj,T) from this result that 
gives the frequency waiting of the functions a. For example, the results of 
Caldeira and Leggett can be re-expressed as 

a^{t -s)= [ —J{uj, T) cosuj{t - s) coth 

J TT 2 

-^J{u!,T)smu!{t- s) (4.67) 

In our case, however, we must first integrate over the angular dependence of k. 
This, however, gives the sum of two terms with Bessel function factors of order 
and 2, themselves dependent on the wavenumber k and the coordinate path 
X(t). In order to define a spectral function, we would have to disentangle the 
t — s and k behaviours, which, with the additional Ji(A;|X(t) — X(s)| factors is 
rather involved. 



4.3.4 Interpreting the imaginary part 

The one magnon coupling treated perturbatively endows the vortex with an 
effective mass and introduces dissipation. In the influence formalism, we expect 
to obtain terms in the effective action of the forward/backward paths inter- 
pretable as particle-like inertial terms. Dissipation arises due to fluctuating 
forces inflicted by scattered magnons on the vortex. We expect the fluctuating 
forces to be accompanied by corresponding damping forces. 

Our one magnon term couples to the vortex velocity and not position as treated 
by Caldeira and Leggett. This is because the vortex is a solution itself of the 
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system, so that all first order variations vanish. The velocity term survives 
because the vortex is to zeroth order a stationary solution. The potential renor- 
malization found earlier going like should here appear as a shift ~ X{t)'^, or 
an inertial term from which we can deduce an effective vortex mass. 

Substituting for /k into the imaginary term yields the phase, including the ad- 
ditional minus sign in (|4.65|l . 

$ = J d^kj\tj\s (X(0e*-^W - Y(t)e^'"Y(*) 

where we define the phase angles via 

F = exp {^{'^-i^)^ 

Performing first the integral over Xk from to tt (refer to an identical calculation 
in the perturbation calculation of multi- vortex mass corrections in section ll:.2.1|l 
to leading order yields, for the term only for conciseness (the other factors 
have the same form), 

^ = -^-^jdkj\tf^ (^X(i) . X(.) Jo(fc|X(0 ~ X(.)l) 

+ X{t)X{s) ((X(t) • eA)(X(s) • ga) (4.70) 

-(X(i) X ga) • (X(s) X ga)) J2(fc|X(0 - X(s)l) j + etc. 

where ga denotes the unit vector connecting X(t) and X(s). 

Integrate by parts in t — s to get two terms, one with two time derivatives in 
X{s) and another with a single time derivative in X{s). Note, we ignore the 
derivatives of the Bessel functions since the extra factor of k makes these higher 
order corrections and we assume the vortex curves slowly to ignore derivatives 
of the unit vectors. The boundary terms from the integration by parts are zero 
for t = s and otherwise unimportant (they don't contribute to the equations of 
motion, being just boundary dependant). Finally, we have 



„ smuj}^{t — s) 
kQ^r^ 

(4.68) 
(4.69) 



Sq^TT 



dk kdtds 



2a2 

X Jo(fc|X(t)-X(s)|) 



sinujoit - s) ^ ccoscjo(t-s) ^ \ 

kQ^ri Q^l ) ■ 

sin.;o(^^-.) ^(^^_ ccosa;og-.) ^(^^^| (4.71) 



X(s) • ( GA(X(t) • ga) - Gi(X(t) • e^))] J2{k\X{t) - X(s)|) + etc. 
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where e_L is a unit vector perpendicular to ga- 

Split this integral into the sine and cosine portions, $ = + $c- Consider first 
the sine integrals. Integrate by parts again in i — s. The non-zero boundary 
terms are 



or equivalently, 



(4.72) 



^Bc^_sfjrrl ^\^^±2^^^_^2^^^^^Rs ^^^^^^ 



lo ^ ' ^Ja^ +rl 

where we've again split the sine integral into $s = + 

These provide an inertial mass term to the effective action of each the forward 
and backward paths. 

The remaining terms, ignoring Bcsscl function derivatives as before. 



J Jo Jo ckQ'^r^ 
Jo{k\X{t) - X{s)\) + etc. 

-^^Jo ''Jo'' ckQ^rl 

Jo(fc|X(t) - X(s)l) + etc. (4.74) 



This is much smaller than the log divergent boundary term. Note that we've 
neglected the smaller still J2 terms. By varying the XX terms with respect 
to X, we would obtain a small third order time derivative term, X, in the 
equations of motion. In an attempt to explain their numerical simulation results, 
Mertens et. aLiii^ii^ artificially introduce a third order term by expanding the 
energy functional assuming both position and velocity as collective coordinates. 
This, of course, is a misapplication of the collective coordinate formalism, where 
a collective coordinate is meant to replace a continuous symmetry that the 
soliton breaks. The freedom they introduced by assuming velocity as a collective 
variable is not actually available in the original problem. 



Consider the cosine term next. We can re-express this damping in terms of 
various damping functions 
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i-T i-t 



J dtj ds(jii{t- s,\X{t) -X{s)\)±(s) -Xit) (4.75) 

+ 7A(i - s, \X{t) - X{s)\)X{s)X{s) ■ eA(X(t) • ga) 

+ 7i(t - s, \X{t) - X{s)\)X{s)±{s) ■ e^{X{t) ■ e^)) + etc. 



where 



2 J Q'^rl 
lA{t - s, A) = / dk k ——- 



2 _ 

7i(i-s,A) =-7A(i-s,A) (4.76) 

Note we cannot perform the k integrals in analytic form due to the loq — ckQ 
argument in the cosine. 

The damping forces depend on the previous motion of the vortex. These memory 
effects appear as averages over Bessel functions — this form is because the vortex 
exists in a 2D system. The first damping term is of the regular form, that is, a 
force acting in the opposite direction to the particle velocity. The next damping, 
7a, is the same as the first if the vortex travels in a straight line, however, for 
a curved path, is dependent on its change in direction. The last damping, 
contributes damping perpendicular to the 7a damping, which, in the case of a 
slowly curving path, is transverse to the vortex motion. 

Comparing with the dissipation results of Slonczewski^^, although we find fre- 
quency dependent dissipation (via the kX{t) coupling in the Bessel functions), 
we do not see any of the same small frequency behaviour predicted by Slon- 
czewski. Likely, his treatment considers a different source of dissipation than 
the contribution considered here. As noted in section 14.2.21 this dissipation 
arises due to the same scattering processes that yield an inertial energy. In 
Slonczewski's treatment, on the other hand, his inertial energy calculation is for 
intermediate distance magnon scattering, while his dissipation arises from far 
field scattering. 



4.3.5 Interpreting the real part 

In the paper of Caldeira and Leggett— , the real part of the influence functional 
is interpreted as the correlation of forces in the classical regime. The real phase 
of their influence functional is 
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r = V coth / dt ds {x{t) - y{t)) cosa;fe(i - s) {x{s) - y{s)) 

'-^ zmuk z Jo Jo 

(4.77) 

which they compare to the contribution of a normally distributed classical fluc- 
tuating force F{t) with correlation {F(t)F{s)) = A{t — s) 

t = ^j\tj^ ds {x{t) - y{t)) A{t - s) {x{s) - y{s)) (4.78) 

Since these terms have the same form, the real part of the influence functional 
must be interpretable as the correlation of forces in the classical regime. 

The real phase of the vortex influence functional is, after substitution for fk, 

r^^Jd'k coth ^ £ dt £ ds (x(t)e*-^(*) - Y{t)e^-^^'^) ■ 

(x(s)e-'-^(^) - Y{s)e-^^-^^^A ■ 

(4.79) 

The integral over (pk can be performed exactly as was done for the imaginary 
part yielding Bessel function pre-factors 

r = !S/d.fccoth^rdJ*d/°^'^^(^--^) 



J 2 Jo Jo kQ^rl 

^X(i) • X(s) Jo(fc|X(t) - X(s)l) + X{t)X{s) ((X(i) • eA)(X(s) • ba) 

-(X(t) X ba) • (X(s) X ba)) J2(A;|X(t) - X(s)|)^ + etc. (4.80) 

Integrating by parts twice to cast this into a similar form to the regular dissi- 
pation term of Caldeira and Leggett 

r = !^ /dA^^oth^ Cdt fds'^^^^^^ 

2a-' J " 2 Jo Jo Q^rl 

^X(t) • X(s)Jo(fc|X(t) - X(s)l) + {(Kit) ■ eA)(X(s) • ba) 

-(X(t) X ba) • (X(s) X ba)) J2(fc|X(i) - X(s)l) ) + etc. (4.81) 
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where, as usual, we neglect derivatives of the Bessel functions since their deriva- 
tives provide higher order corrections. There are additional boundary terms 
depending only on initial and final positions that don't affect the vortex dy- 
namics. 

The real phase can be interpreted as the correlation of forces. However, here 
the fluctuating forces are now vector forces and there are correlations between 
various components of the fluctuating forces. The appearance of the various 
Bessel functions, arising because the vortex is an extended object in 2D, differs 
from the treatment of Caldeira and Leggett because of a different density of 
states of the environmental modes. 



4.4 Discussion of vortex effective dynamics 

Markovian approximation 

We can apply the Markovian approximation as in Castro Neto and Caldeira's 
treatment of solitons"*^. That is, approximate ^{t) — s- 773^ and similarly in the 
force correlation integral, l|4.81|l . giving 



J. MP) 

t M = 



h 

where 



/ dt(X(i) - Y(0) • (X(i) ~ Y(t)) (4.82) 
Jo 



A{P) = S^Jirq^h [ dk k-^ coth (4.83) 
J Q^rf, 2 



Note, all J2 terms disappear in this approximation. 

In this limit, the longitudinal damping coefhcient, (|4.75(l becomes 



(4.84) 



The classical fluctuation-dissipation theorem is now satisfied in the high tem- 
perature limit (coth a: ^ 

i(/3) = 2kBTr] (4.85) 
where T here denotes temperature. 
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This limit corresponds to the hmit where the timescale of interest is much greater 
than the correlation time of the magnons. 



4.4.1 Comparison with radiative dissipation 

The dissipation found in the Markovian approximation can be compared with 
the over-simplified calculation performed using second order perturbation the- 
ory. There, assuming only the emission of a magnon and no inter-magnon 
scattering, we found that the dissipation rate was given by the integral 



where we evaluated this integral in section r4.2.2l Comparing the k dependance 
of this integral with that of the damping coefficient in equation H4.84|l . we find 
they differ only by the (5-function. The dissipation is now dependent on the 
entire magnon spectrum. In the previous calculation, we made the simplifying 
assumption that there were initially no magnons and hence only zero energy 
magnons could be scattered. 

In the Markovian limit, the effective damping force found in the imaginary part 
of the influence functional phase gives roughly the energy dissipation 



Here, the full spectrum of magnons is excited, with probability of finding a 
certain k state weighted by its corresponding Boltzmann factor. Thus, even as- 
suming no vortex inertial energy, we can find scattering between infinitesimally 
spaced k states throughout the spectrum. 



4.4.2 Extending results to many vortices 

The entire treatment can be repeated for a collection of vortices. Assuming the 
vortices are well enough separated to neglect core interactions, the unperturbed 




(4.86) 




(4.87) 
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spin configuration is 



n 



Hot 



i=l 



n 



0, 



tot 



(4.88) 



where denotes the center of the ith vortex. The center coordinates are 
elevated to operators within the collective coordinate formalism. Expanding 
the Lagrangian in terms of this spin field, without magnon interactions, we find 
gyrotropic momentum terms and inter-vortex potentials 



where Ey^i is the unimportant rest energy of the vortices, Pgyro.z = — ^ Xj x 
z is the vector potential giving the gyrotropic force, and the last term accounts 
for inter- vortex interactions. 

The magnon interactions are treated to leading order only — the zero-point en- 
ergy shifts do not affecting dynamics, and higher order dissipation is not treated 
here. There is a one magnon coupling with the vortex velocities X^ to the 
magnons that is integrated over in the influence functional. 

The resulting influence functional now has effective action terms coupling not 
only the forward and backward paths of the same vortex, but also the paths for 
different vortices. Without going through all the details, the general results are 
presented. Were we to neglect all inter-vortex terms, the influence functional 
would simply be the product over each single vortex influence functional. 

Including inter- vortex terms to leading order now, the mass tensor is exactly the 
same as that found using second order perturbation theory (see section f4. 2. l|l 
and, in fact, the calculations here are nearly identical to those. The mass tensor 
is 



There are inter- vortex damping forces behaving essentially like the single vortex 
damping forces: there exist forces longitudinal and transverse to the motion of 




(4.89) 




(4.90) 



where e^, = 



|X.-Xj|- 
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a vortex, however, acting on a second vortex. The damping decreases as a 
function of vortex separation as ~ J^ikrij). This dissipation is thus quite smaU 
when we assume well separated vortices, in keeping with previous calculations 
(refer to the inter- vortex forces calculation in section rrT|l . 

Similarly, in the force correlation integral, we find that the fluctuating forces 
acting on various vortices are inter-correlated. This shouldn't be surprising at 
all: we have damping terms intermingling the motion of vortex pairs so that we 
should therefore expect that the fluctuating forces on these vortices are inter- 
dependent. 

The final effective density matrix propagator becomes 



J(X„Y,;X^,YO = TT / " '2?[X,(i),Y,(t)]expJ(5,[X,(i)]-5,[Y,(t)]) 

exp-^V / / dtdsf VA,,(i-.s)X,(i)-X,(s)Jo(fc|X,(i)~X,(s)|) 
n ^ Jo Jo \ 

+ J2 ^"(i - s)((X,(t) • eAj(X,(s) • eA j 



(X,(t) X gaJ • (X,(s) X gaJ J2(fc|X,(t) - X,(s)|) + etc. (4.91) 



where the force correlations as applied to vortices i and j are 

- = / O^J tolcosukit - s) (4.92) 



The vortex effective action has been redefined to include the inertial mass and 
damping terms 

= I dticl + Y.\^^3^^'^i 

- f ds(Y.l]^{t - s, \X,{t) - X,(,s)|)X,(.) . X,{t) (4.93) 
+ Y.^A{t - s, |X,(t) - X,(s)|)i:,(s)(X,(s) • eA,)(X,(i) • gaO 

i 

+ 7l(^ - s, |X,(t) - X,(s)|)X,(s)(X,(s) • e^.) (X,(i) • e^)) ] 
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where the damping functions arc from (|4.7()|) 



S^Jirq^qj /" , cosu^pft - s) Jo(fcA) 



, dk k „ 
2 J Q 



2 . 

YAt-s,A)=-Y^it-s,A) (4.94) 

Note in the hmit of slow motion and large inter- vortex separation, that is, 
Jo ^ 1 for same vortex terms and all others are negligible, this effective action 
has the same form for each vortex as found in the quantum Brownian motion de- 
scribed by Caldeira and Leggett^, however, with inter-vortex terms introducing 
Coulomb-like forces. 



4.4.3 Frequency dependent motion 

Perhaps a better way of understanding the role of the Bessel function pre- factors 
is to decompose them according to the sum rules 

oo 

J,(A:|x-y|)= J2 Jm{kx)J,+,niky)e'^''+"''>^^--^y^ (4.95) 

7n— — OO 

Denote X^^^ — J„i(fcX^)e*"^'^^X. The effective Lagrangian is transformed to 
^Cl + Y^l dk\ M^iC,„, ■ XL, + i E • Krn+2e'''^' 

rn ij i 

-fds(Y. ^p(* - • Xi,„(t) (4.96) 

i 

where M^^ = jj^gi^a , recalling that the mass tensor can be expressed as an 
integral over k with Bessel function factors (refer to section |4.2.1|I . The new 
damping function is defined as 

= 2 0^2^°"'^°^^''^ 

i /j. \ S'^jTTqf k /, \ 

lAkit - s) = ^ cosa;o(i - s) 

7lfe(t-5)=-7Xfe(i-s,A) (4.97) 
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(4.98) 



(4.99) 

Recall that the density matrix propagator is not simply the product of non- 
interacting forward and backward paths. As in (|4.91|) . we also have damping 
terms coupling the forward and backward paths. 

Thus, we find that the motion of the collection of vortices behaves as interacting 
Brownian particles; however, with frequency dependent damping and fluctuat- 
ing forces. The formalism of Caldeira and Leggett^ can be applied to each 
frequency component, with the added complexity of inter-vortex forces. 

4.4.4 Summary 

A collection of vortices are quantized by considering the small perturbations 
about them. This amounts to including vortex-magnons interactions. We con- 
sidered two couplings in depth: a first order coupling between the vortex ve- 
locity and the magnon spin field, and a second order magnon coupling. All 
vortex-magnon couplings create dissipation via magnon radiative processes. We 
considered only the dissipation due to the first order coupling, first in perturba- 
tion theory and later via the influence functional. Higher order couplings also 
create dissipation, and may, in fact, contribute more significantlyiiS*SS, however, 
these weren't considered here. 

The one magnon coupling creates an inertial energy endowing the vortex with 
an effective mass. A collection of vortices are strongly coupled: in addition 
to the usual inter-vortex forces, there are inter-vortex inertial terms such as 
iMy Xi • Xj that are non-negligible. The zero point energy shift from the two 
magnon coupling is log divergent and, being due to the presence of the vortex, is 
considered the quantized vortex's zero point energy. Note, we did not calculate 
the full effect of this two magnon coupling, only that portion independent of 
magnon populations. This shift was calculated first by considering magnon scat- 
tering in Chapter Inland next in this chapter by simply rewriting the interaction 
in terms of magnon creation/annihilation operators. 



The real part of the influence functional can be re-expressed now as 
^-\T. f '^^ f d^^i - (^^" (^) - Yfc™ it)) ■ (Xi„ {s) - 

+ 5,,A^i{t - s)e'2^-XL„(t) . (^eA,(XI.,„+2(s) • eA.) 

- e_Li • {y^k,m+2{s) • e_Li)^ + etc. in XY and 

where 
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The influence functional reconfirms the effective mass calculations and gives 
explicitly the damping forces and corresponding fluctuating forces responsible 
for dissipation. These act longitudinally and transverse to the vortex motion. 
Again, a collection of vortices are coupled via the damping forces: damping 
forces due to the motion of a first vortex act on a second vortex. Damping 
forces depend on the entire history of the vortex dynamics. 
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Conclusions 



We study the dynamics of a collection of magnetic vortices in an easy plane 
two dimensional insulating ferromagnet. The system is approximated by a con- 
tinuous spin field because we are only interested in the low energy response. 
The vortices interact with magnons via a variety of couplings. The effective 
dynamics bear many similarities to that in the more complex superfluid and 
superconducting vortex bearing systems. 

We reviewed the derivations of the gyrotropic force and the inter-vortex force 

by expanding the vortex action about a stationary superposition of vortex so- 
lutions. We reviewed the inertial mass derivation by calculating vortex profile 
distortions when in motion and showed the equivalence of this method with or- 
dinary perturbation theory. We reviewed magnon phase shift calculations and 
how these phase shifts modify the vortex zero point energy. By rewriting the 
scattering potential in terms of magnon creation and annihilation operators, we 
found an equivalence of the phase shift calculations with the immediate energy 
shift revealed in the second quantized form. 

We suggest a new interpretation of the gyrotopic force as a Lorentz-type force 
with the vortex vorticity behaving like charge, 47reog (in SI units), in an effective 
perpendicular magnetic field, B = -^^piZ, due to the vortex's own out-of-plane 
spins. We rewrite the effective action term giving the gyrotropic force instead 
as a vector potential shift in the vortex momentum. This momentum term was 
then verified by direct integration of the operator generating translations. The 
vector potential possesses gauge freedom, allowed by the same freedom of gauge 
in the Berry's phase. 

We next employed the Feynman- Vernon influence functional formalism, assum- 
ing the vortcx-magnon systems arc initially uncoupled with the magnons in 
thermal equilibrium (thus introducing temperature). The systems interact and 
entangle. The dynamics of the vortices were isolated by tracing over magnons. 
The resulting effective vortex motion is acted upon by longitudinal and trans- 
verse damping forces. Before now, no damping force acting transverse to the 
vortex motion has been suggested in a magnetic system. 

The vortex is a stable solution of the easy plane ferromagnet. As such, when we 
expand about it to quantize magnons in its presence, we find no linear coupling 
between the two fields. However, the vortex is a stationary solution, so that 
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setting it into motion, we find a first order coupling between the magnon field 
and the vortex velocity. 

This lowest order coupling, responsible for endowing the vortex with an effec- 
tive mass, is dissipative and yields effective damping forces acting on a moving 
vortex. The damping forces are accompanied by fluctuating forces that average 
to zero and with time correlations such that the fluctuation-dissipation theorem 
is satisfied in a generalized way. 

We found both longitudinal and transverse damping forces dependent on the 
prior motion of the vortex. For a collection of vortices, the damping forces also 
act between vortices: the motion of a first vortex causes a damping force to act 
on a second vortex. Correspondingly, there are non-zero correlations between 
forces acting on two different vortices. 

The vortex dynamics were described by the propagator of the vortex reduced 
density matrix. The forward and backward paths are coupled, as already de- 
scribed for quantum Brownian motion by Caldeira and Leggett^. The damping 
forces possess memory effects, a common feature in general when describing 
a soliton as a quantum Brownian particle In our two dimensional system, 
however, we found additional Bessel function factors. These considerably com- 
plicate the extraction of a spectral function describing the ensuing Brownian 
motion. By decomposing the vortex motion in a basis of Bessel functions, we 
find that the various frequency components behave as a coupled ensemble of 
quantum Brownian particles. 



5.1 Open questions 

The analogy of a vortex as a charged particle in a magnetic field can be ex- 
tended. For instance, there should be excitations within the gauge field giving 
the gyrotropic momentum. The magnetic field is a result of the out-of-plane 
spins at the vortex center. Perhaps, gauge fluctuations are related to vortex 
core flips. 

Future work on magnetic vortex motion should check the relative importance 
of higher order dissipative couplings. The basic motion of a small collection 
of vortices can be examined now including inertial and damping forces. For 
instance, one could verify the claim of Slonczewski— that damping forces acting 
on a vortex pair only decay circular orbits inward and parallel ones outward. The 
similarities with superfluid vortices should be further examined by attempting 
to calculate the influence functional of a superfiuid vortex and, likewise, the 
Aharanov-Bohm interference effects of magnons passing a moving vortex. 
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Some mechanics 



A classical system is describable by its Lagrangian, which is a function of the 
system coordinates qi and velocities qt 

The action of the system is defined by 



The equations of motion of the system are given by the principle of least action, 
otherwise known as Hamilton's principle, stating that the system evolves from 
initial state {(Zi(0)} to final state {qi{T)} via the path qiit) that extremizes the 
action, S. 

Given that C = C{qi,qi^t), extremizing the action we arrive at the Euler- 
Lagrange equations 



Alternatively, we can describe the system by its Hamiltonian. We transform 
from the Lagrangian to the Hamiltion via a Legrendre transformation 




(A.l) 



d_dC _ d£ 
dt dcii dqi 







(A.2) 




(A.3) 



where we've defined the conjugate momenta Pi defined by 



dC 



(A.4) 




Hamilton's equations are a restatement of (jA.2|l and ljA.4|l 



% _ dH 
dt dpi ' 



dpt 
dt 



dH 

dqi 



(A.5) 
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For example, consider a particle of mass m, position x, residing in a potential 
V{x). The Hamiltonian is simply the total energy of the system 

H{x,p,t)^^ + Vix) (A.6) 

where the conjugate momentum p = mi, as usual. The Lagrangian is then 
given as the difference in kinetic and potential energies 

C{x, i, t) = ^mi^ - V{x) (A.7) 

Application of either the Euler-Lagrange equation or Hamilton's equations yields 
Newton's second law of motion, F = mx, where the force F = —■^V{x). 

Define the Poisson bracket {•, -{q.p for a system with coordinate q and conjugate 
momentum p via 

, , dA dB dAdB 

Note that {q,p}q^p — 1. The above can be easily generalized to a field theory 
by substituting q — + (^{x) and p — > 7r(x) and replacing all simple derivatives by 
functional derivatives. 

Going over to quantum mechanics, to quantize the motion of the system, we 
impose the commutation relations 

[q,p\=ih (A.9) 



In 1925, P. A. M. DiraoiS observed that proper quantum mechanical relations 
followed under the substitution 

in 

In a spin system, using coordinate (p a-nd conjugate momentum Scos9 we can 
verify directly the classical version of [Si, Sj] = ihetjkSk, that is, 

{Si, Sj}cf,,scose — SijkSk (A. 10) 



where we define S = ^(sin^ cos 0, sin sin 0, cos 0). However, spin being an 
essentially quantum concept, we must bear in mind that when speaking of spin 
directions given by {(j),9), we mean the spin state of highest probability to be 
found in that direction. 



A.l. Imaginary time path integral 
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A.l Imaginary time path integral 

Consider a system in thermal equilibrium at temperature t. If we decompose 
the system Hamiltonian into a set of eigenstates £,n{x) with cigenenergies En, 
then the probability of observing the system in eigenstate n is proportional to 
e ''B^ where ks is Boltzmann's constanti^. The density matrix for this system 
is 

p{x',x) = Y,U^')Cr,{x)e-^''- (A.ll) 

n 

where f3 — [ksT)^^ . 

Compare this with the quantum propagator decomposed into this same basis: 



iHT 

K(x',T;x,0) =(x'|exp 

ri 

= ^(.T'|6.)exp-^(C„|x) (A.12) 



h 

iE„T 



We see that the density matrix is formally identical to the propagator corre- 
sponding to an imaginary time interval T = —iph. In fact, if we consider the 
equation of motion of the density matrix found by taking the derivative of (jA.lip 
with respect to 

^Y.^nUx')Cn{^)e-^''- (A.13) 

Recall that En£,n{x') — H^n{x'). If we understand to act only on x', we 
can write 

-^^^^n.'p{x',x) (A.14) 

We know how to evaluate the propagator as a path integral for simple Hamilto- 
nians involving only the system coordinates and their conjugate momenta. For 
example, for the Hamiltonian 

H = - — —^ + V{x) (A.15) 
2m ax^ 



the solution over an infinitesimal time period e is 
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' ^' = V ^ n [J—T^ '''^[-^)) ^^-^'^ 

which can be verified by direct substitution into 
_^ dKix'^T;x,0) ^ 

Now, under an infinitesimal interval in the density matrix ie the solution is given 
by e — + — ie 



= e/H) = ^^exp-- (y^-^ + j j (A.18) 

which can be verified by direct substitution into (jA.14|l . 

Stringing many of these solutions together for successive intervals of time ac- 
cording to 

p(.',./,', = /..'V(.'..",^'W."..;ffl (A.19) 

for intermediate x" at (3, we obtain a path integral formulation of the den- 
sity matrix which is simply an imaginary time version of the propagator path 
integral, that is, with the substitution T — > —iph. 
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Appendix B 

Quantization of classical 
solutions 

Suppose we have a particle described by position x residing in a potential V{x). 
Classically, the particle follows a path x{t) that satisfies Newton's second law of 
motion. In quantum theory, the particle is no longer described by its position 
X, but by its wavefunction '4^{x) giving a probability distribution of finding the 
particle at position x. If the energy is conserved, the wavefunction can be 
decomposed into energy eigenstates, ilinix), obeying Schroedinger's equation 

Hil;^ = EnMx) (B.l) 

where H is the Hamiltonian of the system, quantized by elevating the position 
and momentum variables to operators. 

As preparation for a description of the quantization of a soliton solution, con- 
sider some of the finer points of quantization of classical particle solutions. 

For the potential shown in Figure IB. II there are three extrema and hence three 
stationary classical solutions. The absolute minimum, a; = a is the classical 
ground state, having the lowest attainable energy. 

In quantum mechanics, according to the uncertainty principle, a solution is not 
allowed to have zero momentum and a fixed position. Thus, even in its ground 




J I I I 

a b c X 



Figure B.l: An illustrative potential of a one dimensional particle. 
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state there are fluctuations. Expanding V{x) in a Taylor scries, to lowest order 
the potential is harmonic about the minimum and we have simple harmonic 
excitations with frequency = V"{x = a) and energies 

Er, = V{a) + {n+^)nuj (B.2) 

The ground state energy becomes Eq = V{a) + \hu:. The additional ^hu) is the 
zero-point energy due to quantum fluctuations. 

The solution a; = c is a second stable classical solution. Quantum mechanically, 
there are again fluctuations about this solution that give a similar excitation 
spectrum. In this case, however, since this is an excited state, there are possibly 
tunneling processes that relax the state to its ground state about x = a. In a 
field theory, this stable excited state is the analogue of a soliton solution (while 
the tunneling processes are analogous to instantons). However, for the mag- 
netic solitonic solutions considered here, these excited states belong to separate 
topological sectors of the solution space so that there is effectively an infinite 
energy barrier to the ground state. 

The classical solution x = 6 is unstable and would thus correspond to an imag- 
inary frequency. There are hence no set of quantum levels formed about it. 

Another interesting analogy to consider is the case of a constant potential, 
V{x) = V . In that case, there is no clear choice of minimum about which to 
expand and, should we attempt to, we would find everywhere w = 0. Of course, 
in quantum mechanics, the proper solutions to consider are the plane waves e'*^^ 
with energies 

En = V + ^{nKf (B.3) 

where ?ifc„ = p„ are the momenta of these states. In field theory, we find a zero 
frequency mode, or Goldstone mode, for every broken continuous symmetry. 
Further, for each of these broken symmetries, we find a corresponding conserved 
momentum, analogous to the conserved p„ in the particle case. 



B.l Quantizing soliton solutions 

In field theory, quantizing a soliton follows analogously to the regular quanti- 
zation of a classical solution. The language is changed somewhat however. For 
instance, the ground state of the particle, a; = a, is quantized to a hierarchy 
of simple harmonic excitations. In field theory, we call the absolute potential 
minimum the ground state, or the vacuum state. The hierarchy of perturbative 
excitations are interpreted as mesons or quasiparticles. In our system, these are 
the magnons. 



B.l. Quantizing soliton solutions 
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When wc expand about the solitonic excited state (analogous to the second 
minimum, x = c), generally the quasiparticles are modified by the soliton pres- 
ence. In the simple particle case, this corresponds to the general case where 
V"{a) V"{c). In the particle case, the hierarchy of simple harmonic states 
are interpreted as excited states about the minima. 

In a field theory, the quasiparticles are generally extended states and, in the 
presence of a soliton, are shifted but still extended. In some cases, the soliton can 
trap a few quasiparticlc modes. These bound modes arc interpreted as soliton 
excited states. The remaining, extended states are interpreted as unshifted 
quasiparticles, while all energy shifts due to the soliton are attributed to the 
zero-point energy of the quantized soliton. 

The soliton acts perturbatively on the extended states, it itself being localized 
in space, as a scattering center. Asymptotically far from the soliton center, the 
quasiparticles are simply phase shifted. Suppose that the relative phase shift 
between the incoming and outgoing waves is 6{k), a function of the wavevector 
k. 

By enforcing periodic boundary conditions^ on both the unperturbed wavevec- 
tor k and the scattered wavevector q 



Lkn =2n7r 

Lqn - 5{qn) =1n-n (B.4) 

we fix the allowed k and q values. In the L oo limit, these allowed values 
merge to a continuum and the sum over fc-states is replaced by an integral 



L 
2^ 



dk 

oo 



The energy correction to the soliton solution, taken as the modification to the 
zero point energy of the vacuum, is thus, noting that = w(k ^ j) , 



AE=hY^io{q)-Lu{k) 



1 r ^ 



tOr alternatively, we could enforce fixed boundary conditions forcing k to be 7r-periodic 
rather than 27r-periodic. In that case, the asymptotic waveform must be modified from a 
plane wave to a cosine wavefunction and we find that the phase shift is also changed by a 
factor of 2. Thus, either set of boundary conditions is equivalent. 
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found by expanding to first order in S. 

In addition to small corrections to the quasiparticle continuum, the soliton might 
bind discrete levels in the quasiparticle spectrum. Those with a; = are due to 
a continuous symmetry broken by the soliton solution. These modes are dealt 
with using collective coordinates. There can also be w ^ discrete modes. 
These are interpreted as soliton excited states. For an example of these, see the 
quantization of the quantum kink of the (f)"^ theoryi^S — the magnetic vortex does 
not have any such excited states. 

B.1.1 In a path integral formalism 

Using path integrals, we can find the excitation spectrum of a system by taking 
the trace of the system's quantum propagator. We first review the simple case 
of a regular particle in an external potential and then generalize to field theory. 

Semiclassical approximation for a single particle 

The propagator of a single particle starting in position qa at time and ending 
in position qb at time T is 



where H{q,p) — ip^ 
Next, we set qa = qb 



G{T) 



where {4>n} denote a complete orthonormal set of eigenstates of H. This yields 
an expression giving the excitation spectrum of the Hamiltonian. 

For a particle in a potential V{x) with a minimum at a; = a;o, the classical 
solution is simply qd = xo- Expanding the potential in a power series about 
this solution 



K{qb,T;qa,0) =< qb\e-'"^^^\qa > (B.6) 

- V{x), and, for simplicity, we've set m = 1. 

= go and integrate over the endpoint of the periodic orbit 



dqo<qo\e-'"^/^\qo> 
- fdqoY,< > e-^^"^/'' < 0„|go > (B.7) 

n 



V{x) = V{qci) + V'{q,i){q - qa) + ^V" {qa){q - qdf + OiAx^) (B.8) 



B.l. Quantizing soliton solutions 
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the second term is zero since qd is a minimum of V{x). The action expanded 
about this classical solution, q{t) q^ + q'{t), is now 

S[qit)]^-Vixo) + dt^q'^ ~]^uj\'^ (B.9) 

where uP = V"{xq), assumed positive (i.e. the classical solution is stable). 
Note, at this point, the boundary conditions of the periodic path are still not 
generally satisfied so that the new perturbed solution must now satisfy q'(0) = 
q'{T) =qo- qd- 

The semiclassical approximation amounts to neglecting the 0{Ax^) and higher 
order terms. But the terms in q' are just the action of a simple harmonic 
oscillator. To evaluate the path integral 

GsHoiT) = JdqoJ I?[g'(t)]e^^[''Wl (B.IO) 

we expand again about the simple harmonic oscillator classical solution satis- 
fying the appropriate boundary conditions. You may ask why we didn't im- 
mediately go from the beginning action and expand in a Taylor series and 
approximate there. Although that would have proceeded identically, the addi- 
tional step helps clarify what to do when expanding in a field theory admitting 
classical soliton solutions. The classical solution is now 

q'^i = A cos Lut + B sin cut (B.ll) 
where the boundary conditions give 

A =qo - qd 

A cos uiT + B sin ojT ^qo - qd (B.12) 



Evaluating this second classical contribution to the action gives 

S[q'J^~2uiqo-qdr'^^^^ (B.13) 

The complete path integral becomes, setting y{t) = q'(t) — q'^i and noting y[t) 
now has the boundary conditions y{Q) = y{T) = 0, 

G(T) = y dgoe-^^("")-^""('''-«^"'''?^ j V[y{t)]e^ 'i^vi-l^-^')v 

1 r^. 



Solving for the determinant of the remaining action — ^ /p dty{-^ + uj'^)y we 



solve the relevant Jacobi equation (^ + uj'^)y — with initial conditions 
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y{0) = and y'{0) = 1. This gives the prefactor 



2Trih sin wT 
Evaluating the qo integral, the final result is 



G(r) ^ 



(B.14) 



2i smLoT/2 



1 - e-*'^^ ^ ' 



n=0 



giving the excitation spectrum = ?ia;(n + 5) as expected. 



Semiclassical approximation in field theory 

This follows almost identically to the single particle case, with just a few tech- 
nical points needing clarification. Suppose we have a field theory in 1+1 dimen- 
sions with the Lagrangian density 

C,{x,t) = \{d^4>f-U{,l^\ (B.16) 



Assume [x) is a stationary extremum of this system. Expanding the action 
about this solution, (j) ^ cp' + (j)ci 

S = S,i + \jdxj dtid.d^'f - ^^^</.'' (B.17) 

Next, we integrate by parts to replace {d^(f)')'^ <i>'{—^ + ■^)4>' ■ 

Assuming now that (l)'{x,t) is separable, i.e. (j)'{x,t) = f{x)g{t), we solve for 
the eigenvalues of the spatial portion 

-^ + ^^)m-)=-.V.(.) (B.18) 



Assume that the fr{x) eigenfunctions form an orthonormal basis. Expressing 
the general solution (f>'{x,t) = fr{x)gr{t) so that the integration measure 
becomes Hr ■^[5'''(*)]' action becomes 



B.l. Quantizing soliton solutions 
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Iff 

^=2 1^^ '^^^^'■^^^3r{t)^{--Q^-0jl,)fr'{x)gr'{t) 

r t' 

= Y.\\ dtgrm--^ - uj'Mt) (B.19) 

by the orthonormality of the fr{x). Thus, the problem has separated into a 
product on r of equivalent single particle problems 

G{T) = e*^"' Yl (^J D[5^(t)]e*-^*''^(*)(-^-'"')'?-(*)^ (B.20) 

which wc know how to solve from the previous section. 
The only remaining manipulation required is to note that 

J|^e-''^'^'-("'-+5) = e-^S'-^'^'-("'-+5) (B.21) 

r rir {rir} 

where {ur} denotes a set of integers rir- 
B.l. 2 Collective coordinates 

Suppose the soliton exists in a system with translational symmetry. The soli- 
ton itself is a localized entity, and hence breaks this symmetry. The soliton 
must choose arbitrarily what coordinate to center on. This is an example of 
spontaneously broken symmetry. 

This symmetry introduces to the quasiparticle spectrum a zero frequency mode 

associated with the soliton. While to first order presenting no problems, should 
we continue in the perturbative expansion, the energy denominators would de- 
velop artificial singularities. 

In perturbing about the soliton solution, rather than as done previously via 

oo 

(f) =(f)0 + ^ an{t)%l)n{x) 
n=0 

7 1 oo 

=<^o + ao(0^ + $^an(t)V'n(a;) (B.22) 

n— 1 

where the n = mode is the translation mode, rewrite the expansion as 

oo 

<^ = <^o(a;-X(i)) + ^a„(i)V'n(a;) (B.23) 
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where X{t) is the coUective coordinate associated to the translational invariance. 
This is completely equivalent if we expand 4'o{x — X{t)) to first order in X{t) 
and identify ao{t) with —X{t). 

Rewriting the Lagrangian in terms of this expansion, the potential terms, being 
translationally invariant by assumption, does not depend on X{t). The kinetic 
term depends only on X{t). 

We can introduce conjugate momenta to X{t) and to the an{t), denote these P 
and 7r„, and transform to the classical Hamiltonian 

C30 

H = PX{t) + ^ 7r„a„(t) - L (B.24) 

n=l 

Quantizing the soliton now follows exactly as quantization of a regular particle: 
we impose commutation relations on the various degrees of freedom 

[X, P] =ih 

[a„,7r„] ^ih (B.25) 

The quantized quasiparticles have a zero-point energy shifted by ^ ^hdcun that 
is attributed instead to the quantized soliton. That is, if the vacuum quasipar- 
ticle zero-point energy is ^ 5?ii-^n, while in the presence of a soliton becomes 
^ ■^Ti{ujn + Sujn), the soliton is said to have the zero-point energy ^ ^hSujn 
while the quasiparticles are considered unchanged—. 
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Spin path integrals 



Consider a spin system with Hamiltonian H. The propagator for this spin 
to evolve from state Vli at time t = to il/ at time i = T is 

i^(%,T;r!„0) - (r!/|exp-Ji7T|a) (C.l) 

Inserting — 1 resohitions of the identityi^ 

2s + 1 



/ dVL\n)(n\ = 1 (C.2) 

47r J 

where a lower case s denotes the dimensionless spin (whereas, S = ?is), gives 

(l]7V-i|---|f^i)(r!i|exp-Jif^|r!o) 

n N 



where fc = denotes the initial state and k = N the final state. Define e = T/N. 
Expand the exponential 

Keeping terms to linear order in e, the H term can be approximated at equal 
times: define H{tk) = (f2fc+i |ff |f2fc). Re-exponentiate the bracketed term to 
exp~j-eH{tk). 

The overlap of two coherent states, flk and flk+i is^ 

{nk+^\nk) = (l±Ih±ii^^ (C.3) 

where 
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^ = 2 tan- ftan ( ^^±1^) cos ^0.^, + e,)\ _ 

and where ^ is a gauge dependent phase that we can ignore. The pre-factor is 
1 to leading order and the phase can be approximated such that 

(flfe+i |rifc) = exp ^-zse^^iy— ^ cos 6k^ (C.5) 
AU together, letting N — + oo, we find the spin path integral 



K{nf,T;n„0) = J V[n{t)]exp (-is J dt^{t) cos 6 (t) ~ H{t)\ (C.6) 



Note that there are no spurious boundary terms as there are in the stereographic 
representation using z and z*, as found, for example, by Solari^-. 



C.l The semiclassical approximation 

Evaluation of spin path integral is non-trivial as evidenced by the series of papers 
suggesting various corrections. Klauder— discussed the spin path integral in 
terms of conjugate variables (j) and 5 cos 6* and first addressed the semiclassical 
approximation applied to the spin path integral. He claimed that to evaluate 
properly the trace of the propagator obtaining the excitation spectrum, real 
valued periodic orbits are required. However, simple counting of degrees of 
freedom, given two equations of motion (one for (f) and another for Scos9) with 
two initial and two final conditions, results in an overdetermined system. In 
fact, we are also trying to simultaneously specify both x and p at each boundary, 
disallowed by the familiar uncertainty principle. 

Kuratsuji and Mizobuchi^^ note this overdeterminacy and claim only one of 
{xi,Xf} or {pi,Pf} needs specifying, the other being fixed by the equations of 
motion. 

Solari— finds an additional pre-factor 

exp^ [ dtA{t) (C.7) 
2 Jo 

where A(t) is a time-dependent operator appearing in the action zA{t)z* where 
z is the spin coherent state in the stereographic projection. We won't worry 
about this correction since in our treatment there is no such term in the action. 
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Various authors have even claimed that the spin path integral can only be 
properly evaluated discretely. However, a continuous version is reliable with the 
proper additional phase of Solari, as argued by Stone et. al.^. 

Below, we review the usage of the Jacobi equation for evaluating the path in- 
tegral of a regular particle, then generalizing to the path integral over a field. 
Finally, we derive the analogous Jacobi equation for a spin path integral, fol- 
lowing closely the work of Kuratsuji— . 



C.1.1 Coherent state path integral 

In the classical limit, a spin coherent state \^{t)) can be interpreted simply as 
a spin lying along the direction 17 (i). 

The spin path integral, including the trace over periodic orbits, can be written 



as^ 



Jn(o)=n{T) 

This step is analogous to the quantum perturbations about a soliton solution. 
We now are solving for the quantum propagator for these perturbations. 

Let (j)ci{t) and Oci{t) be a classical solution of this action (analogous to the simple 
harmonic oscillator solutions of the single particle case). Attempting to impose 
periodic boundary conditions results, in general, in an over-determined system 
of equations. Instead, we set only (/)c/(0) = <j)ci{T) allowing the equations of 
motion to fix boundary conditions for 0ci{t). 

Expanding cj) = (j)ci{t) +x{t) and SO — SOci{t) + y{t), the action becomes to sec- 
ond order variations (neglecting higher orders in keeping with the semiclassical 
approximation) 

S = Sa- j dt (^xy + i {A{t)x^ + 2B{t)xy + C{t)y^)^ (C.9) 

where A{t) = 0, B{t) ^ and dt) = ^^i^. 

In the discrete version^, introducing the small timestep e, we complete the square 
in yk to obtain 



tPor the moment considering a single spin - the generalization to a field of spins follows 
identically to the treatment in Appendix 1^ 

tin arriving at this expression, note that in the discrete version there is actually an average 
of y(t) yk+Vk-i under careful analysis gives boundary terms as found by Solari»S. 

We neglect these terms and approximate y^.. 
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' k=l 



( Bk Axk- Xk-\)\ 



Notice wc use the extra integration over the periodic orbit coordinate On to 
integrate over all N j/fe's; whereas, we only use N — 1 integrations over the x^s. 
Impose the boundary conditions xq = xn = 0- 

For the general case, where we do not have the additional integration over 
boundary conditions, we must introduce this additional integration as an aver- 
aging over the final coordinate. This doesn't change the physics since this final 
coordinate is necessarily fixed by the equations of motion anyway. 

The N Gaussian integrals over yk give the pre-factors Ylk=i iM^- "^^^ com- 
plete expression becomes 



G(T) = lim — 
^ ^ jv^oo V 27r 



jV 



N 



.fe=l 



Kk=l 



exp 



fe=l 



^k 



-exlak (C.U) 



where ak=^Ak-^ + i{§)^. 



The problem becomes that of solving for the determinant of the {N — 1) x (A'' — 1) 
matrix 



—I 



V 



where = -tak + + ic^- 

Re-express the product of pre-factors from the yk Gaussian integrals as 

det . ieCi 



( 



N 

n 

fe=i 



2-K 

ieCk 



C.l. The semiclassical approximation 
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and noting that det{AB) = det{A) det{B), we multiply the two matrices to yield 



/ Caeai -1 



iCiedet 



V 



CN-lCaN-2 —1 



Denote the determinant of the submatrix ending in the fc'th row and column by 
Dk- We can then write down the recursion relation 



-Dfe = eCk+iakDk-i 77^£>fe-2 



2 d_ fB 
dt [c 



Dk-i 



Letting D/. be a function of ke, this can be rewritten 



Dk — 2-D/c-i + -Dfc-2 _ (Cfc+i — Ck){Pk-\ — -Dfc-2) 



or in a continuum limit 



df^ ~ C dt dt V C ^ dt\C 



(C.12) 



The initial conditions on D can be found directly from the first and second 
submatrix determinants 



Di = iC\eC2ea\ 

D2 — Di ( ~ ^ ~ C3 _ 

iCi ( C3ea2C2eai - — — 62eai 
(^2 



giving, in the limit of e ^ 0, -D(O) = and -D(O) = iC{0). But this is equivalent 
to the system of equations from the original formulation 



I = A. + By 
dt 



dx 
Itt 



-Bx - Cy 



(CIS) 
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with initial conditions a;(0) = and y(0) = --1 after eliminating y{t) and setting 
ix{t) — D{t). Thus the required determinant is ix{T). 



C.1.2 Spectrum of a ferromagnetic plane of spins 

The action of a ferromagnetic plane of spins with easy plane anisotropy in a 
continuum limit is 

s^s j j dt(j)9 - - (^-0v20 - e\/^e + -2 j (C.14) 

where the various constants are as defined in Chapter |2| 

Choose a set of spatial eigenfunctions such that — > —k^ and J ^fk'fk = 
S^k'-k). 

Thus, the integration measure becomes a product over k states, now decoupled, 
leaving within the time integral of the action (note S was factored out into the 
integration measure) 



where to = ckQ. The periodic classical solutions can be written 

till )=a( T^"' , ) (C.16) 

with the periodicity condition on 0fe(i) imposing identical conditions on A and 
B as in (|B.12|I . with qo^qd 4'kQ- Note that the periodicity condition was pre- 
viously (j){x, 0) — 4>{x, T) — (jjQ] however, after the transformation to diagonalize 
the equations in fc, each coefficient must now be periodic and integrated 
over. 

The classical action for these periodic orbits becomes 

S. = -2^0L^^^^ (C.17) 

The perturbed action has the same form as the original linearized action above, 
ljC.14p . Calling the small perturbations in (j), x, and those in 6, y, we need a 
solution such that a;(0) — and y{0) — 1 (the change of sign here arises from 
linearizing cos 9 —>■ —9 in the Berry phase term). This corresponds to 



y{t) I \ coswfct 



(C.18) 
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and the determinant evaluates to ix{T) = sinwfcT. 
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